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A  TECHNIQUE  TO  MEASURE  THE  VOLUME  OF  ELASTOMERS  AS  A 
FUNCTION  OF  TEMPERATURE  AND  PRESSURE  WITH  AN  ACOUSTIC  PYCNOMETER 


I.  INTRODUCTION 


This  report  describes  the  use  of  an  acoustic  pycnometer  to  measure 
the  volume  of  a  spherical  elastomer  sample  as  a  function  of  temperature 
and  pressure.  The  measurement  of  the  volume  of  elastomer  samples  is 
needed  in  determining  the  bulk  modulus  of  elastomers. 


The  bulk  modulus  of  a  material  is  defined  as  the  ratio  of  a  tensile 
or  compressive  stress,  trlaxial  and  equal  in  all  directions  (e.g. , 
hydrostatic  pressure),  [11  to  the  relative  change  in  volume  it  produces. 
For  a  spherical  sample  this  trlaxial  compressive  stress  is  depicted  as 
a  squeezing  of  the  spherical  sample  by  hydrostatic  pressure,  in  which  a 
change  in  volume  occurs  without  a  change  in  shape  (See  Fig.  1)  [2],  The 
static  bulk  modulus  is  defined  as 


B(P,T) 


AP 

AV/Vg(P,T)  ’ 


(1) 


where 

B(P,T)  *  the  static  bulk  modulus  at  pressure  P  and  temperature  T  . 

AV  ■  the  change  in  volume  of  the  sphere 

V  (P,T)  *  the  volume  of  the  spherical  sample  at  P  and  T. 
s 

To  calculate  the  static  or  isothermal  bulk  modulus  using  Eq.  (1), 
a  precise  method  to  measure  the  volume  of  the  sample  V  (P,T)  and 
correspondingly  the  change  in  volume  of  the  sample  AV  is  needed. 
Manuscript  submitted  August  14,  1981. 
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The  dynamic  bulk  modulus  of  a  material  can  be  defined  as  the  ratio 
of  a  harmonically  varying  acoustic  pressure  to  the  corresponding  change 
in  volume  it  produces,  multiplied  by  the  volume  of  the  sample  of  the 
material  at  the  appropriate  static  pressure  P  and  temperature  T.  The 
expression  for  the  dynamic  bulk  modulus  of  the  sphere  in  Fig.  1  is 


K(P,T,f) 


p(f  ) 

Av/Vg(P,T)  ’ 


(2) 


where 

K(P,T,f)  =>  the  dynamic  bulk  modulus  of  the  sphere  at  static  pressure 
P,  temperature  T,  and  frequency  f 
p  *  the  varying  uniform  acoustic  pressure  with  frequency  f 
iv  »  the  change  in  volume  of  the  sphere  due  to  the  acoustic 
pressure  p. 

Persons  who  work  in  the  field  of  hydroacoustics  are  interested  in 
the  dynamic  bulk  modulus  of  elastomers.  Elastomers  are  used  in 
transducers  to  isolate  the  electrical  components  from  seawater. 

Another  application  of  elastomers  in  hydroacoustics,  is  in  their  use 
as  hull  coatings  on  submarines.  Therefore,  a  precise  experimental 
method  to  determine  the  dynamic  bulk  modulus  is  of  importance.  One 
way  to  obtain  an  expression  for  the  dynamic  hulk  modulus  of  an  elastomer 
sample  is  from  the  theory  describing  the  scattering  of  sound  bv  an  elastic 
sphere.  in  brief,  if  one  uses  the  theorv  ror  Kavleteh  scattering  as  the 
basis  for  lvoamlc  bulk-modulus  determination,  the  following  arguments 

apply. 
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An  expression  containing  the  dynamic  bulk  modulus  of  an  elastomer 
is  obtained  by  solving  the  boundary- value  problem  for  the  scattering  of 


sound  by  an  elastic  sphere  [3]  composed  of  the  elastomer  that  is  immersed 
in  a  fluid.  If  one  looks  at  the  scattering  of  sound  in  a  confined  chamber 
where  a  volume-expander  sound  source  (shown  in  Fig.  2)  produces  a  harmonic, 
uniform  acoustic  pressure,  then  the  expression  for  the  scattered  pressure 
is 


.3  2 

PiVe 

3r 


(3) 


where 

P  *  the  scattered  acoustic  pressure 
s 

p^  *  the  incident  acoustic  pressure 

k^  =  2Trf/c^,  the  wave  number  in  the  fluid 

c^  ”  the  speed  of  sound  in  the  fluid 

ag  *  the  radius  of  the  spherical  elastomer 

”  the  dynamic  bulk  modulus  of  the  fluid 

■  the  dynamic  bulk  modulus  of  the  spherical  elastomer 

r  ■  the  distance  from  the  center  of  the  spherical  elastomer 
to  the  receiver. 

Equation  (3)  is  true  for  nearfield  (k^r  <  1)  Rayleigh  scattering 

(k,a  <<:  1)  of  the  zero  mode.  The  zero  mode  is  the  first  term  in  the 
j  e 

series  expansion  of  the  solution  to  the  boundary-value  problem  and 
physically  represents  a  uniform  acoustic  pressure  field  or  breathing 


SPHERICAL 


ACOUSTIC 

SOURCE 


ELASTOMER  SURROUNDING 
SAMPLE  FLUID 


RECEIVER 


RIGID 

CHAMBER 


Fig.  2  —  Scattering  of  sound  by  a  spherical  elastomer  sample 
in  a  confined  chamber 
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mode.  The  breathing  mode  is  the  mode  of  vibration  that  corresponds 
to  a  harmonically  varying  hydrostatic  pressure.  If  the  chamber  is 
sufficiently  small,  only  the  zero  mode  is  present. 


If  an  acoustic  scattering  experiment  is  performed  as  shown  in 
Fig.  2,  the  scattered  pressure  is  not  measured  directly.  The  total 
pressure  is  measured  by  the  receiver  and  is  given  by  the  expression 


ptotal  "  Pi  +  Ps’ 

where 


(4) 


ptotal  *  t*'e  tota^  Pressure. 


An  expression  for  the  ratio  of  K^/Kg  in  the  confined  chamber  in 
Fig.  2  can  be  obtained  if  three  measurements  are  made:  a  measurement 
with  the  chamber  filled  only  with  fluid,  a  measurement  with  a  spherical 
elastomer  sample  present  in  the  fluid-filled  chamber,  and  a  measurement 
with  a  rigid  scatterer  present  in  the  fluid-filled  chamber.  By  measuring 
the  total  pressure  Ptota2  with  no  sample  present  in  the  chamber,  the 
incident  pressure  expressed  by  Eq.  (4)  can  be  found.  In  any  nearfield 
experiment,  the  distance  r  from  point  of  observation  to  the  sample  must 
be  taken  into  account.  This  is  accomplished  by  introducing  the  expression 
for  the  scattered  pressure  for  a  rigid  sphere  when  placed  in  the  same 
position  in  the  confined  chamber  in  Fig.  2.  The  expression  [4]  for  the 
scattered  pressure  for  a  rigid  sphere  is 
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where 


p  =  the  scattered  pressure  from  a  rigid  sphere 
r ,  s 

a^  =  the  radius  of  the  rigid  sphere. 

Using  Eqs.  (3),  (4),  and  (5),  one  can  obtain  expressions  for  the 
total  acoustic  pressures  for  scattering  by  a  spherical  elastomer  sample 
and  by  a  rigid  sphere  in  the  confined  chamber  in  Fig.  2.  The  results 
are 


t,e 


,  2  3 

k3ae  Kf 

1 - ■=— £  (i - L 

3r  U  Ke  * 


(6a) 


and 


P 


t,r 


.2  3 
k3ar 
3r 


(6b) 


where 

p  *  the  total  pressure  measured  when  a  spherical  elastomer 
t ,  6 

sample  is  present  in  the  confined  chamber 
^  =  the  total  pressure  measured  when  a  rigid  sphere  is  present 
in  the  confined  chamber. 

The  ratio  of  the  bulk  modulus  of  the  fluid  to  the  bulk  modulus  of 
the  spherical  elastomer  sample  is  obtained  by  combining  the  results 
given  in  Eqs.  (3),  (4),  and  (6),  assuming  that  p^  is  the  same: 


where 

V^_  =  the  volume  of  the  rigid  sphere 

=  the  volume  of  the  spherical  elastomer  sample. 

Equation  (7)  expresses  the  ratio  of  the  dynamic  bulk  modulus  of  the 
fluid  in  the  confined  chamber  (refer  to  Fig.  2)  to  the  dynamic  bulk 
modulus  of  the  spherical  elastomer  in  the  chamber.  To  calculate  the 
ratio  K^/K^,  using  Eq.  (7),  as  a  function  of  temperature  and  pressure, 
the  volume  of  the  spherical  elastomer  sample  must  he  known  as  a 
function  of  temperature  and  pressure.  A  confined  chamber  is  analyzed 
because  one  is  interested  in  determining  the  ratio  K^/K^  as  a  function 
of  pressure  and  temperature.  The  chamber  shown  in  Fig.  2  permits  one 
to  simulate  the  effect  of  ocean  pressures. 

The  volume  of  a  spherical  elastomer  sample  can  be  measured 
accurately  by  means  of  an  acoustic  pycnometer.  The  background  for  the 
acoustic  pycnometer  used  evolved  from  the  work  of  Corsaro,  Jarzynski, 
and  Davis  [5],  These  investigators  describe  the  use  of  an  acoustic 
densitometer  to  measure  the  change  in  volume  of  polyethylene  oxide  as 
a  function  of  temperature  and  pressure  with  a  precision  of  3  parts  in 
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where 

V  =  the  volume  of  the  rigid  sphere 

Vg  =  the  volume  of  the  spherical  elastomer  sample. 

Equation  (7)  expresses  the  ratio  of  the  dynamic  bulk  modulus  of  the 

fluid  in  the  confined  chamber  (refer  to  Fig.  2)  to  the  dynamic  bulk 

modulus  of  the  spherical  elastomer  in  the  chamber.  To  calculate  the 

ratio  Kf/Ke»  using  Eq.  (7),  as  a  function  of  temperature  and  pressure, 

the  volume  of  the  spherical  elastomer  sample  must  be  known  as  a 

function  of  temperature  and  pressure.  A  confined  chamber  is  analyzed 

because  one  Is  Interested  In  determining  the  ratio  K,/K  as  a  function 

r  e 

of  pressure  and  temperature.  The  chamber  shown  in  Fig.  2  permits  one 
to  simulate  the  effect  of  ocean  pressures. 

The  volume  of  a  spherical  elastomer  sample  can  be  measured 
accurately  by  means  of  an  acoustic  pycnometer.  The  background  for  the 
acoustic  pycnometer  used  evolved  from  the  work  of  Corsaro,  Jarzynski, 
and  Davis  [5].  These  investigators  describe  the  use  of  an  acoustic 
densitometer  to  measure  the  change  in  volume  of  polyethylene  oxide  as 
a  function  of  temperature  and  pressure  with  a  precision  of  3  parts  in 
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II.  ACOUSTIC  PYCNOMETER  THEORY 


The  acoustic  pycnometer*  is  an  instrument  which  is  used  to  precisely 
measure  the  volume  of  an  elastomer  sample.  In  this  section  the  theory 
used  in  the  acoustic  pycnometer  measurement  is  outlined.  First, 
the  principle  of  operation  of  the  pycnometer  is  described.  Second,  the 
pertinent  equations  are  derived. 


The  acoustic  pycnometer,  shown  in  Fig.  3a,  is  composed  of  a  two- 
part  Invar  container,  held  together  with  screws  and  sealed  with  neoprene 
O-rings.  Within  the  container  there  is  a  sample  chamber  connected  to 
a  bore  by  a  narrow  canal.  At  the  base  of  the  bore  is  an  ultrasonic 
quartz  transducer.  The  transducer  is  a  resonant  device,  which  is  1.905 
cm  in  diameter,  0.75  cm  in  thickness  and  operates  at  a  frequency  of 
4  MHz.  The  entire  pycnometer  is  filled  with  mercury  (shown  in  black 
in  Fig.  3a)  under  vacuum.  The  sound  path,  in  the  bore  of  the  pycnometer, 
is  terminated  by  a  resonant  float  to  flatten  the  mercury  meniscus.  The 
float  design  will  be  discussed  later. 


*The  acoustic  densitometer  of  Corsaro,  Jarzynski,  and  Davis  was  obtained 
from  the  Naval  Research  Laboratory  in  Washington,  DC,  but  because  of  its 
existing  condition  was  unusable  when  this  research  was  undertaken. 

Several  modifications  had  to  be  made  in  order  to  use  it.  Because  of 
these  extensive  modifications  and  because  the  modified  device  was  used 
in  a  different  way  than  by  Corsaro,  Jarzynski,  and  Davis,  its  name  was 
changed  to  the  acoustic  pycnometer  to  distinguish  it  from  the  original 
densitometer.  The  new  name  is  appropriate  because  a  pycnometer  is  an 
instrument  used  to  measure  the  volume  or  the  density  of  a  material.  On 
the  other  hand,  the  word  "densitometer"  might  imply  a  connection  with  the 
kind  of  instrument  used  in  photography  when  there  is  none. 
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Fig.  3a  —  Acoustic  pycnometer  with  a  spherical  sample  present 


An  electronically  generated  4-MHz  pulse  is  used  (see  Fig.  4)  to 
excite  the  transducer  at  the  base  of  the  bore  of  the  pycnometer  causing 
the  transducer  to  transmit  an  ultrasonic  4-MHz  pulse  through  the  column 
of  mercury  in  the  bore.  This  pulse  is  reflected  by  the  float  and 
travels  back  through  the  column  of  mercury  where  the  reflected  pulse  is 
received  by  the  transducer  and  is  also  reflected  by  the  transducer  back 
into  the  mercury-filled  column.  The  process  repeats,  producing  a  sequence 
of  echoes,  until  an  electronically  controlled  repetition  time  has  elapsed 
and  another  pulse  is  generated  by  the  electronic  equipment.  This  electronic 
pulse  excites  the  transducer,  and  another  transmitted  pulse  is  produced 
at  the  transducer.  Again  a  series  of  pulse  echoes  occurs.  This 
series  of  pulse  echoes  is  commonly  called  an  echo  train.  Figure  5  is  a 
photograph  taken  of  an  oscilloscope  display  of  an  echo  train. 

When  a  sample  is  placed  in  the  sample  chamber  of  the  pycnometer 
(see  Fig.  3a),  the  change  in  volume  of  the  sample  caused  by  any  change 
in  temperature  and  pressure  produces  a  change  in  the  volume  of  the 
mercury  in  the  bore.  The  height  of  the  mercury  column  is  determined  as 
follows.  The  time  it  takes  for  a  4-MHz  sound  pulse  to  travel  from  the 
transducer  in  the  base  of  the  bore  through  the  mercury  column  to  the 
float  and  return  to  the  transducer  is  measured.  This  time  is  called 
the  time-of-f light.  This  time-of-f light  is  determined  by  measuring  the 
difference  in  time  between  two  successive  pulse  echoes  on  an  oscilloscope 
(as  shown  in  Fig.  4).  By  multiplying  the  time-of-f light  by  the  speed 
of  sound  in  mercury,  the  length  of  the  mercury  column  for  any  temperature 
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Fig.  4—4  MHz  driving  pulse,  timebase  =  2  ps 
Pulse  amplitude  =  600  volts,  peak-to-peak 


Fig.  5  —  Pulse-echo  train,  timebase  =  50  /js  ampli¬ 
tude  of  the  first  echo  in  the  train  is  approximately 
=  0.8  volts  peak-to-peak. 


and  pressure  can  be  calculated  by 


L  "  2  'lig  td  »  (8) 

where 

L  *  the  length  of  the  mercury  column 
c^  *  the  speed  of  sound  in  mercury 
t^  =  the  difference  in  time  between  two  successive  pulse-echoes. 

The  measurement  of  the  time-of-f light  will  be  discussed  in  detail 
in  Appendix  A.  The  volume  of  the  mercury  in  the  bore  of  the  pycnometer 
is  found  by  multiplying  the  length  given  by  Eq.  (8)  by  the  cross- 
sectional  area  of  the  bore.  The  measured  volume  V  is  corrected 
to  account  for  the  change  in  dimension  of  the  bore  produced  by  changes 
in  pressure  and  temperature.  The  volume  of  the  mercury  in  the  bore  is 
given  by 

VBORE  =  2  cHg  Cd  irrb  exPl2^“vT  “  BVP^3^’  ^ 


where 

r,  *  the  radius  of  the  bore 
b 

Oy  a  the  volume  coefficient  of  thermal  expansion  of  the 
pycnometer  material 

*  the  volume  coefficient  of  isothermal  compressibility  of 
the  pycnometer  material. 
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The  units  of  the  temperature  and  pressure  are  respectively  in  degrees  C 
and  MPa  and  the  pressure  used  is  gauge  pressure.  These  units  for 
temperature  and  pressure  will  be  used  throughout  this  paper.  Equation 
(9)  is  derived  in  Appendix  C. 

Mercury  is  used  as  a  working  medium  in  the  acoustic  pycnometer 
because  the  speed  of  sound  and  compressibility  of  this  liquid  are  well 
documented.  The  equation  for  the  speed  of  sound  in  mercury  as  a  function 
of  temperature  and  pressure  is  [6] 

c^F.T)  -  1460.0  -  0.460T  +  0.210P.  (10) 

The  units  of  the  sound  speed  in  Eq.  (10)  are  m/s.  The  volume  of  an 
elastomer  sample  can  be  calculated  if  the  volume  of  mercury  in  the  bore 
is  determined  in  two  cases:  a  measurement  is  made  with  a  sample  present 
in  the  sample  chamber  of  the  pycnometer  (shown  in  Fig.  3a)  and  a  second 
measurement  with  no  sample  present  in  the  sample  chamber  (see  Fig.  3b). 
The  measurements  taken  with  no  sample  present  in  the  pycnometer  are 
referred  to  as  a  calibration  run.  A  calibration  method  is  used  to 
account  for  the  compressibility  of  the  pycnometer,  which  is  a  significant 
effect  in  the  experiment. 

To  obtain  an  expression  for  the  volume  of  a  sample  when  using  the 
acoustic  pycnometer,  two  cases  must  be  analyzed:  the  calibration  run 
and  the  sample  run.  Writing  down  the  equations  for  the  volumes  of 
mercury  in  the  pycnometer,  in  the  configuration  shown  respectively 
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in  Figs.  3a  and  3b,  one  obtains 


V  ip  T}  »  V  +  v  +V  -V 

HgV  ’  '  CHAM  CANAL  BORE  SAMPLE 


(11a) 


and 


V  fP  Tl  =*  V  +V  +V 

Hgvr,w  VCHAM  CANAL  BORE  * 


(lib) 


where 


V  (P,T)  =  the  volume  of  mercury  in  the  pycnometer  during  the 
h8 

sample  run 

I 

VHg(P,T)  *  the  volume  of  mercury  in  the  pycnometer  during  the 
calibration  run 


v 

CHAM 

=  the 

volume 

of 

the  chamber 

VCANAL 

■  the 

volume 

of 

mercury  in  the 

canal 

V 

BORE 

=  the 

volume 

of 

mercury  in  the 

bore 

during 

the 

sample  run 

V' 

BORE 

=  the 

volume 

of 

mercury  in  the 

bore 

during 

the 

calibration  run 

v  SAMPLE 

=  the 

volume 

of 

the  sample. 

To  obtain  an  expression  for  the  volume  of  the  sample,  one  subtracts 
Eq.  (lib)  from  Eq.  (11a).  The  result  is 


/  m  “V  J.  111  »■ 

'sample  “  BORE  ”  BORE  Pu  (T.) 

ng  r 


(12) 


where 


m 


■  the  mass  of  mercury  used  in  the  pycnometer  during 


the  calibration  run 


m  =  the  mass  of  mercury  used  in  the  pycnometer  during  the  sample 


run 

PHg(T^)  ■  the  density  of  mercury  at  the  filling  temperature  and  the 
filling  pressure  P^. 

In  Eq.  (12)  the  difference  V  '  -  V„  in  the  volumes  of  mercury  in  the 

Hg  Hg 

pycnometer,  which  arose  from  the  subtraction  of  Eq.  (lib)  from  Eq.  (11a), 
has  been  replaced  by  the  difference  in  the  masses  of  mercury  divided  by 

f 

the  density  of  mercury.  The  difference  in  V,,,,-..  and  VDnn„  at  a  given 
temperature  T  and  pressure  P  is  due  to  the  change  in  volume  of  the 
sample  and  to  filling  the  pycnometer  with  mercury  to  different  heights 
during  the  sample  run  and  during  the  calibration  run.  Therefore, 


vbore(p,t) 


V„  _  + 
HgB 


AV 


SAMPLE 


(13a 


where 


V  (T  ) 
HgB'  f ' 


AV 

SAMPLE 


■  the  volume  of  the  mercury  in  the  bore  at  temperature  Tf 
and  pressure  P^ 

■  the  change  in  volume  of  the  sample. 


VgoRE  is  equal  to  VB0RE  plus  a  volume  of  mercury  due  to  filling  the 
bore  of  the  pycnometer  to  a  different  height: 


VunD*(P,T)  “  VU„D  +  V (P,T)  , 


(13b 


V(P,T)  =  the  additional  volume  of  mercury  in  the  bore  due  to 

filling  the  pycnometer  to  different  heights  during  the 
sample  run  and  during  the  calibration  run. 

I 

The  difference  in  VB0RE  and  Vg0RE  is  obtained  by  subtracting  Eq.  (13a) 
from  Eq.  (13b).  The  result  is 


V„^(P,T)  -  V  (P ,T ) 


BORE' 


BORE' 

1  SAMPLE' 


(14a) 


■  (P.T)  -  VSAMpLECPo,To)  -  V(P,T), 


where 

VSAMPLE^Po,To^  =  the  volume  of  the  sample  when  it  is  placed 

into  the  pycnometer  (as  measured  by  Archimedes' 
principle)  at  temperature  Tq  and  ambient  pressure 

P  . 

o 

It  should  be  noted  that  the  pressures  P^  and  Pq  both  refer  to  ambient 
atmospheric  pressure.  Also  the  temperatures  Tq  and  Tf  were  both 
approximately  25°C.  By  rearranging  Eq.  (14a),  one  obtains  an  expression 
for  the  volume  of  a  sample  as  a  function  of  pressure  and  temperature: 

VSAMPLE(P,T)  "  VSAMPLE(Po,To)  , 

(141 

+  vbore<p't)  -  vbore<p-t)  +  V(P'T)- 

Before  Eq.  ( 1 4b)  can  be  used,  an  expression  is  needed  for  V(P,T) 
in  terms  of  known  or  measured  quantities.  This  expression  is  found  in 
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the  following  way.  At  the  pressure  Pq  and  the  temperature  Tq,  Eqs.  (12) 
and  (14b)  when  subtracted  give 


V‘W 


Hg 


m  ~  111  —  v  (v  T  ) 

(P  ,T  )  SAMPLE^  o'1  o'** 


(15) 


Equation  (15)  expresses  the  difference  between  the  volumes  of  the 
mercury  contained  in  the  bore  of  the  pycnometer  during  the  calibration 
run  and  during  the  sample  run  at  temperature  and  pressure  Pq.  As 
mentioned  earlier,  the  sample  volume  V  (P  ,T  )  in  Eq.  (15) 

bAflr  Lh  O  O 

is  determined  by  Archimedes'  principle  at  ambient  temperature  Tq  and 
pressure  Pq,  before  being  placed  into  the  acoustic  pycnometer.  By  using 
Eq.  (15)  and  the  equation  of  Grindley  and  Lind  [7],  which  gives  the  volume 
of  mercury  as  a  function  of  temperature  and  pressure,  one  obtains  for 
V(P,T)  the  expression 


V(P,T) 


D(P.T) 

»<Po,To) 


-  V  (P  T  ) 

SAMPLE'-  o’  o ' 


(16) 


where 

D(P,T)  =  1  +  dxT  -  d2P  -  d3TP  +  d^P2  +  d5TP2 

dj  =■  1.821  ^  10“A 

d2  -  3.153  x  io“5 
d3  »  5.19  x  io~8 
dA  -  4.54  x  io"9 
d5  -  1.19  x  io"11. 
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From  the  previous  development  of  the  theory,  it  can  be  seen  that 
the  volume  of  a  sample  can  be  determined  as  a  function  of  temperature 
and  pressure  by  Eqs.  (14b)  and  (16).  By  using  the  calibration  method, 
the  compressibility  of  the  pycnometer  can  be  taken  into  account. 

Previous  investigators  attributed  the  compressibility  of  the  pycnometer 
to  be  due  solely  to  the  compressibility  of  the  Invar  container  material 
acting  as  a  monolithic  block  of  material.  The  calibration  method  requires 
accurate  measurement  of  the  initial  sample  volume  prior  to  being  inserted 
into  the  sample  chamber  of  the  pycnometer.  This  technique  also  requires 
accurate  determination  of  the  amounts  of  mercury  used  in  the  pycnometer 
during  the  calibration  run  and  during  the  sample  run.  However,  all  of 
the  measurements  necessary  for  the  success  of  the  calibration  method 
are.  made  by  standard  laboratory  procedures. 

III.  ACOUSTIC  PYCNOMETER  EXPERIMENTAL  METHOD 

The  acoustic  pycnometer  can  be  used  to  determine  the  volume  of 
elastomers  as  a  function  of  temperature  and  pressure  after  two  preliminary 
steps  have  been  completed.  First,  the  spherical  elastomer  samples  need 
to  be  fabricated.  Second,  one  needs  to  obtain  the  initial  volume  of 
the  samples  at  ambient  pressure  Pq  and  a  known  temperature  Tq 

VSAMPLE(Po,To)' 

The  pycnometer  is  assembled,  and  the  transducer  is  put  into  place. 

The  assembled  pycnometer  is  then  placed  in  a  vacuum  chamber  where  it  is 
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to  be  filled  with  mercury  (see  Fig.  6).  The  mercury  filling  is  done  by 
using  a  burette  from  which  mercury  is  drained  into  the  sample  chamber  of 
the  pycnometer  while  the  pycnometer  is  kept  under  vacuum.  Previous  in¬ 
vestigators  at  our  laboratory  determined  the  volume  of  mercury  drained 
into  the  pycnometer  by  reading  the  volume  of  mercury  used  from  the 
graduated  scale  on  the  burette.  In  the  more  accurate  method,  used  in 
this  investigation,  the  mass  of  the  mercury  used  in  filling  the  pycnometer 
was  determined  by  weighing  and  the  tabular  value  of  the  density  of 
mercury  was  used  to  calculate  the  volume  of  mercury  [8],  The  volume  of 
the  mercury  placed  into  the  pycnometer  is  calculated  by 


mi  "mf 

vv 


(17) 


where, 

m^  =  the  mass  of  mercury  contained  in  the  burette  prior  to  filling 
the  pycnometer 

mj  =  the  mass  of  mercury  contained  in  the  burette  after  the 
pycnometer  has  been  filled. 

The  mass  in  Eq.  (17)  is  obtained  by  weighing  an  initial  amount  of 
mercury,  which  is  then  poured  into  the  burette  shown  in  Fig.  6.  Once 
the  pycnometer  is  completely  filled,  the  remaining  mass  m^.  of  mercurv 
in  the  burette  is  reweighed  and  the  volume  of  the  mercury  in  the  pycnometer 
is  calculated  by  Eq.  (17).  Since  mercury  vapor  is  extremely  toxic, 

exercised  when  performing  the.  mercury  filling. 


extreme  caution  was 


23 


A  major  source  of  experimental  error  can  occur  when  filling  the 

acoustic  pycnometer  with  mercury.  Any  mercury  lost  upon  filling  will 

result  in  an  error  when  determining  the  amount  of  mercury  in  the  acoustic 

3 

pycnometer.  For  instance,  a  0.2-cm  loss  of  mercury  when  filling  the 
pycnometer,  in  the  calibration  run  or  sample  run,  would  result  in  an 
error  of  1.9  percent  in  determining  the  amount  of  mercury  in  the  acoustic 
pycnometer.  This  error  is  significantly  higher  than  the  error  experienced 
when  the  measurements  were  made.  A  second  problem  occurs  when  filling 
the  pycnometer  with  mercury.  Erttrapped  air  is  a  major  concern  when  using 
the  acoustic  pycnometer.  Since  the  sample  chamber  is  cylindrical  in 
shape,  if  the  pycnometer  cell  is  not  completely  level  when  filling, 
air  could  be  entrapped  in  the  top  of  the  sample  chamber.  To  lessen  the 
possibility  of  the  entrapment  of  air,  a  conical  sample  chamber  would 
seem  to  be  a  more  logical  design.  If  it  has  a  conical  shape,  even  if 
the  sample  chamber  of  the  pycnometer  is  not  level  when  the  device 
is  filled,  air  will  not  be  trapped  as  easily. 

After  filling,  the  pycnometer  is  then  ready  to  be  interfaced  with 
the  electronic  equipment  by  way  of  a  waterproof-connector  assembly 
called  the  electronics-interface  module.  This  is  shown  in  Fig.  7.  The 
interface  module  and  a  teflon  ring  make  the  seal  for  the  pressure 
vessel.  The  threaded  pressure-vessel  top  acts  as  a  backing  for  the 
interface  module  to  compress  the  teflon  ring  so  as  to  seal  the  pressure 
vessel.  The  interface  module  also  acts  as  a  feed-through  to  bring 
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Fig.  7  —  Pycnometer  pressure  and  temperature  control  system 


the  electric  wires  from  the  pycnometer.  Inside  the  pressure  vessel, 
out  to  the  electronic  equipment.  The  pressure  vessel  Is  placed  Into 
a  large  water  bath  for  nominal  temperature  control.  Precise  temperature 
control  Is  obtained  by  connecting  a  circulating  bath  to  the  copper  colls 
surrounding  the  pressure  vessel  to  maintain  the  pycnometer  at  the 
temperature  desired.  A  pressure  pump  allows  one  to  control  the 
static  pressure  P  to  which  the  elastomer  sample  Is  subjected. 

Pycnometer  measurements  are  made  as  a  function  of  pressure  at 
constant  temperature.  When  a  desired  temperature  T  is  obtained  at 
ambient  pressure  Pj,  the  volume  of  mercury  In  the  bore  Is  determined  by 
the  time-of-f light  measurement  described  earlier.  The  pressure  is  now 
increased  to  some  new  value  P 2  using  the  pump,  and  the  temperature  is 
allowed  to  return  to  T.  After  the  temperature  has  stabilized,  another 
time-of-f light  measurement  is  made.  This  process  Is  repeated  at  selected 
pressure  values  to  70  MPa,  and  then  this  process  Is  reversed  and  measure¬ 
ments  are  made  at  the  same  pressure  values  back  to  ambient  pressure. 

This  constitutes  a  pressure  cycle  at  temperature  T.  The  selected  pres¬ 
sures  at  which  elastomer  sample  volume  measurements  were  made  were  as 
follows:  0,  1,  2,  4,  6,  8,  10,  15,  20,  25,  30,  40,  50,  60,  and  70  MPa. 
Finer  pressure  increments  were  used  at  low  pressures  (0  to  10  MPa) 
because  the  low-pressure  behavior  of  the  pycnometer  was  desired.  Several 
of  these  pressure  cycles  at  constant  temperature  were  done  for  each  of 
the  volume  measurement  made  during  the  calibration  and  sample  runs. 
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The  data  obtained  from  the  acoustic-pycnometer  measurement  are 
time-of-fllght  as  a  function  of  pressure  at  constant  temperature.  This 
tlme-of-f light  is  the  time  difference  between  two  successive  pulse 
echoes,  shown  In  Fig.  5,  and  Is  equal  to  the  round-trip  time  of  flight 
in  the  bore.  This  time  difference  is  measured  by  electronic  means  as  a 
frequency  which  Is  the  reciprocal  of  the  time  of  flight  of  the  sound 
pulse  In  the  mercury  column  in  the  bore.  A  more  detailed  description  of 
the  electronic  equipment  Is  given  in  Appendix  A. 

The  data  from  the  acoustic-pycnometer  measurements  are  processed  in 
the  sequence  shown  in  Fig.  8.  The  pressure,  temperature,  and  frequency 
data  are  entered  into  a  Hewlett-Packard  Model  9825  Calculator,  which  is 
used  to  compute  the  volumes  of  the  mercury  in  the  bore  during  the  cali¬ 
bration  and  sample  runs,  using  Eq.  (9).  These  results  are  entered  into 
a  curve-fitting  computer  program  RAGT.FTN  on  a  PDP-11/45  computer,  which 
is  used  to  obtain  an  expression  for  the  volume  of  the  mercury  in  the 
bore  as  a  function  of  pressure  at  constant  temperature.  The  computer 
program  RAGT.FTN  uses  the  method  of  least-squares  to  fit  an  empirical 
expression  to  the  data  for  the  volume  of  the  mercury  in  the  bore,  during 
the  calibration  run  and  the  sample  run.  This  empirical  expression  is 
needed  to  perform  the  sample-volume  calculation  using  Eq.  (14b).  The 
computer  program  PYCNOM.FTN  uses  the  results  of  the  program,  RAGT.FTN, 
to  calculate  the  volume  of  the  sample  according  to  Eq .  (1 4b).  The 
-exults  of  these  an.. lyses  are  given  in  Section  V.  The  programs  RACT  'T'' 
■mii  ire  written  in  FORTRAN  and  make  use  of  the  PDP-11  /4S 


INPUT  P,T,F  FROM 
PYCNOMETER  INTO 
HP-9825 


Fig.  8  —  Acoustic-pycnometer  data-processing  flowchart 
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digital  computer  at  the  Naval  Research  Laboratory's  Underwater  Sound 
Reference  Detachment  (NRL-USRD).  The  computer  programs  used  are  given 
in  Appendices  C,  D,  and  E. 

IV.  ENGINEERING  CONSIDERATIONS 

Volume  measurements  initially  were  attempted  with  an  existing 
pycnometer  previously  tried  by  several  experimenters.  The  acoustic 
pycnometer  used  in  this  investigation  required  several  modi f leaf  ions 
before  consistent  experimental  results  were  obtained.  The  major 
modifications  were  the  construction  of  a  thermistor  plug,  a  transducer- 
crystal  mount,  and  a  resonant  reflector  float.  These  modifications  were 
made  as  the  need  for  each  arose  during  the  course  of  the  experiment. 

Figure  9  shows  the  acoustic  densitomer  of  Corsaro,  Jarzynski,  and 
Davis  [5].  These  Investigators  encountered  difficulties  when  making 
measurements  at  pressures  below  7  MPa.  They  attributed  their  difficul¬ 
ties  at  low  pressure  to  the  presence  of  air  in  the  sample  chamber  of  the 
densitometer,  which  entered  when  the  device  was  filled  with  mercury  even 
though  filling  was  carried  out  with  the  device  in  a  vacuum  chamber. 

Also  in  their  design,  the  transducer  in  the  base  of  the  bore  was  backed 
with  a  spring  to  hold  it  tightly  against  the  bottom  of  the  bore  so  as  to 
prevent  mercury  from  leaking  around  the  quartz  crystal  and  shorting  it 
out.  Moreover,  the  float  they  used  on  top  of  the  mercury  meniscus  in 
the  bore  was  a  cylinder  with  a  conical  top.  The  idea  behind  using  the 
conical  top  on  the  float  was  to  provide  a  tapered  impedance  so  that 
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Fig.  9  —  Acoustic  densitometer 


little,  or  no,  reflection  of  the  sound  pulse  would  occur  from  the  top  of 
the  float  back  into  the  mercury  column.  These  investigators  used  a 
platinum  resistance  thermometer  for  monitoring  the  temperature  in  the 
pressure  vessel.  The  theory  Corsaro,  Jarzynski,  and  Davis  developed  for 
the  acoustic  densitometer  system  took  into  account  only  the  compressi¬ 
bility  of  the  Invar  container,  as  if  it  were  a  solid  metal  block,  and  of 
the  mercury  within  the  densitometer.  No  calibration  of  the  system  was 
performed  to  measure  the  behavior  of  the  densitometer  with  the  sample 
absent.  A  volume  determination  of  a  polyethylene-oxide  sample  was  done 
by  performing  only  a  sample  run.  This  neglected  the  possibility  of 
experimentally  determining  the  behavior  of  the  densitometer  as  a  function 
of  temperature  and  pressure.  In  view  of  the  significance  of  the  effect 
determined  in  the  experiment  reported  here,  omitting  the  calibration 
procedure  was  a  serious  shortcoming  in  the  technique  of  the  earlier 
Investigators. 

Previous  investigators  at  this  laboratory  replaced  the  platinum 
resistance  thermometer  used  by  Corsaro,  Jarzynski,  and  Davis  by  two 
thermistors.  One  thermistor  was  located  on  the  outside  of  the  densito¬ 
meter  close  to  the  Invar  container.  The  second  thermistor  was  placed 
inside  the  sample  chamber  of  the  densitometer,  and  the  electrical  leads 
were  brought  outside  the  sample  chamber  through  a  glass-to-metal  seal. 
This  glass-to-metal  seal  was  screwed  into  the  top  of  the  densitomer 
directly  above  the  sample  chamber.  The  addition  of  the  second  thermistor 
assured  that  thermal  equilibrium  was  attained  when  the  two  thermistors 


were  indicating  the  same  value  of  temperature.  When  the  existing 
acoustic  densitometer  was  initially  assembled  and  filled  with  mercurv, 
it  was  observed  that  mercury  was  leaking  around  the  screw  thread*  of  the 

glass-to-met  al  sea’.  'rh<  s  1o'>k  made  the  existing  densitom^r  . ''•’bln. 

A  thermistor  plug  was  made  to  replace  the  glass-to-metn!  sen1  arrangement. 
This  had  a  two-fold  effect:  thermal  equilibrium  could  V  secured  because 
two  thermistors  were  still  used,  and  this  thermistor  plug  provided  a 
better  way  to  fill  the  pycnometer  with  mercury.  Previously,  the  densito¬ 
meter  was  filled  by  draining  mercury  from  a  burette  down  into  the  bore 
where  it  landed  on  the  quartz  crystal.  Sometimes  this  operation  shattered 
the  quartz  transducer  in  the  base  of  the  bore.  Filling  in  this  way  also 
permitted  air  to  be  trapped  in  the  upper  corners  of  the  sample  chamber. 

By  pouring  the  mercury  through  the  thermistor-plug  hole,  the  sample 
chamber  could  be  directly  and  completely  filled  with  mercury.  The 
thermistor  plug  Is  shown  in  Appendix  G. 

To  ensure  a  better  contact  between  the  electrode  and  the  quartz 
transducer  crystal,  an  electrode  plug  was  designed  (see  Appendix  G)  to 
mount  the  electrode  and  transducer  at  the  base  of  the  bore.  This  re¬ 
placed  the  spring-backed  transducer  arrangement  in  the  original  acoustic 
densitometer.  This  transducer  electrode  plug  was  made  of  a  machinable 
ceramic  to  hold  the  transducer  electrode  in  place  without  shorting  it 
out  against  the  Invar  container.  An  O-ring  was  used  to  make  a  seal 
between  the  transducer  crystal  and  the  electrode  plug.  Use  of  this 
electrode  plug  resulted  In  echoes  of  higher  amplitude,  a  well-defined 
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baseline  on  the  received  signal  with  better  quality  pulse  echoes,  and 
more  pulse  echoes  present  in  the  echo  train.  The  electrode  plug  .also 
decreased  the  possibility  of  mercury  flowing  around  the  crystal  and 
shorting  out  the  electrode  because  of  the  O-ring  seal  used. 

The  float  shown  in  Fig.  9  was  replaced  by  the  resonant  float 
shown  in  Fig.  3a.  When  using  the  original  conically  topped  float  and 
observing  the  pulse  echoes  on  an  oscilloscope,  one  found  that  the  base¬ 
line  of  the  signal  was  not  well  defined,  the  amplitudes  of  the  pulse 
eclioes  were  not  very  large,  and  there  were  few  echoes  in  the  echo  train. 
It  appeared  that  this  was  due  in  part  to  an  interference  caused  by  the 
signals  propagating  through  the  float  material  and  reflecting  off  the 
upper  portion  of  the  float.  In  the  new  design,  this  back-face  reflection 
was  used  to  advantage  to  enhance  the  received  signal.  A  resonant  float 
was  made  in  the  shape  of  a  hollow  cylinder,  with  the  top  open  and  the 
lower  end  having  a  thickness  of  one  half  the  wavelength  of  the  sound 
wave.  In  the  material  used  for  the  float,  at  the  driving  frequency  of 
4  MHz.  If  the  reflector  Is  a  half-wavelength  thick,  it  permits  the 
reflected  signal  from  the  inner  surface  of  the  cylindrical-float  base  to 
add  in  phase  with  the  reflection  from  the  base  of  the  float.  This 
Increases  the  received  signal  amplitude  and  consequently  makes  the 
received  signal  more  distinguishable  from  the  baseline.  The  thickness 
necessary  for  the  base  of  the  float  can  be  determined  from  the  operating 
frequency  of  the  transducer  in  the  bore  and  the  speed  of  sound  in  the 
material  used  for  the  float.  The  thickness  is 
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where 


d  =  |c  /f 
2  m 


(18) 


d  =  the  thickness  desired  for  the  base  of  the  float 

c  =  the  speed  of  sound  in  the  float  material 
m 

f  =  the  operating  frequency  of  the  transducer. 

The  float  is  made  of  stainless  steel  for  which  the  bulk  sound  speed 
c^  =  5790  m/s.  As  stated  earlier,  the  operating  frequency  is  4  MHz. 
Therefore,  the  desired  thickness  should  be  0.072  cm  or  0.028  in.  It  was 
found  that  this  float-base  thickness  caused  the  reflected  signal  from 
the  inner  surface  of  the  cylinder  base  to  add  up  in  phase  with  the 
reflected  signal  from  the  base  of  the  float.  This  gave  a  well-defined 
baseline,  an  increased  signal  amplitude,  and  a  greater  number  of  pulse- 
echoes. 

Additional  problems  arose  with  the  existing  densitometer  due  to 
the  nature  of  the  operating  frequency  and  the  high  voltages  encountered 
during  the  experiment.  Electrical  cross  talk  was  abundant  when  the 
system  was  first  used.  Initially,  a  strain-gauge  bridge  with  a  pressure 
transducer  was  used  to  measure  the  static  pressure  in  the  system.  The 
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This  cross-talk  signal  also  appeared  on  the  thermistor  leads,  which  in 
turn  fed  the  signal  into  the  electronic  digital  multimeters 
used  to  monitor  the  temperature  by  measuring  the  resistance  of  the 
thermistor.  As  a  consequence,  the  multimeters  were  also  being  burned 
out.  To  solve  this  problem,  several  capacitors  were  placed  in  parallel 
across  the  inputs  to  the  digital  multimeters  to  short  out  the  cross¬ 
talk  signal  before  it  entered  the  multimeters.  These  shorting  capacitors 
would  have  no  effect  on  the  direct  current  that  flowed  through  the 
thermistor.  In  addition,  the  pulse  amplitude  of  the  electronic  system 
that  generated  the  4-MHz  acoustic  pulse  was  turned  off  when  temperature 
readings  were  being  made.  Similarly,  when  the  pulse  amplitude  was 
turned  on  during  the  time-of-flight  measurement,  the  multimeters  were 
turned  off. 

Once  these  modifications  were  made,  the  acoustic  pycnometer  could 
be  used  to  determine  the  volume  of  elastomer  samples. 

V.  EXPERIMENTAL  RESULTS 

The  acoustic  pycnometer  was  used  to  measure  volume  of  two  spherical 

elastomer  samples  approximately  2.5  cm  in  diameter.  The  two  elastomer 

samples  used  were  a  sample  of  butyl-252  and  a  sample  of  type-W  neoprene. 

The  volumes  of  these  elastomer  samples  were  measured  by  Archimedes' 

principle  at  atmospheric  pressure  and  25°C  before  being  placed  into  the 

sample  chamber  of  the  acoustic  pycnometer.  The  volume  of  the  sample 

3 

of  butyl-252  was  8.8433  cm  ,  and  the  volume  of  the  sample  of  type-W 
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neoprene  was  8.7774  cm  at  25°C.  A  discussion  of  Archimedes'  principle 
is  given  in  Appendix  E. 

As  stated  earlier,  the  three  measurements  taken  when  using  the 
acoustic  pycnometer  are  frequency,  pressure,  and  temperature.  Frequency 
is  the  primary  measurement  necessary  to  obtain  the  time-of-flight 
through  the  mercury  in  the  bore  of  the  pycnometer.  Pressure  and 

temperature  are  auxiliary  measurements.  The  frequency  measurement 

4 

can  be  made  to  two  parts  in  10  .  The  precision  in  the  pressure  measure¬ 
ment  and  the  temperature  measurement  is  respectively  ±0.1  MPa  and  ±0.2°C. 
A  random-error  analysis  is  given  in  Appendix  B.  The  experimental 
results  are  calculated  by  computer  and  plotted  on  a  Tektronix  Model  4662 
flat-bed  plotter. 

The  acoustic  pycnometer  was  calibrated  at  two  temperatures  (10°c 
and  25°C)  and  as  a  function  of  pressure  (to  70  MPa).  The  pycnometer 
calibration  runs  are  presented  in  Figs.  10  and  11.  The  points  shown  In 

I 

the  figures  are  the  measured  values  for  V  the  volume  of  the  mercury 

in  the  bore,  calculated  using  the  HP-9825  program.  The  curves  passing 
through  the  data  points  were  determined  by  the  curve-fitting  routine, 
RAGT.FTN.  Each  had  the  form  of  a  sixth-degree  polynomial.  The 
computer  program  RAGT.FTN  determines  an  expression  for  the  volume 
of  the  mercury  in  the  bore  as  a  function  of  pressure  at  constant 
temperature  by  the  method  of  least  squares.  Bata  were  taken  on  two 
pressure  cycles  at  10°C  and  25°C.  All  of  the  calibration  data  for  a 
given  temperature  ere  analyzed  together,  and  then  the  sixth-degree 
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polynomial  was  passed  through  the  data.  At  pressures  below  2  MPa,  in¬ 
consistent  results  were  obtained.  This  was  apparent  from  the  lack  of 
precision  observed  in  the  data  when  measurements  were  made  at  low 
pressure.  Also,  the  curve-fitting  routine  could  not  fit  the  measured 
values  below  2  MPa.  Therefore,  all  measured  values  below  2  MPa  were 
discarded. 

When  the  measurements  were  made  in  the  case  of  the  sample  runs,  the 
data  For  the  volume  of  the  mercury  in  the  hore  were  analyzed  cycle  by 
cycle,  when  fitting  the  measured  values  to  a  sixth-degree  polynomial. 
Again,  two  pressure  cycles  were  taken  at  temperatures  of  10°C  and  25°C; 
but  instead  of  lumping  these  data  together  and  fitting  a  sixth-degree 
polynomial  to  all  the  data,  a  sixth-degree  polynomial  initially  was  fit 
through  the  volume-  vs  -pressure  data  for  each  cycle.  This  cycle-by¬ 
cycle  analysis  was  done  to  investigate  the  possibility  of  hysteresis  in 
the  elastomer  samples.  The  term  hysteresis,  in  this  sense,  implies  the 
sample  would  not  return  to  its  original  volume  after  pressure  cycling. 
This  hysteresis  was  not  observed  in  the  sample  of  butyl-252.  The  type-W 
neoprene  sample  showed  some  slight  evidence  of  the  hysteresis  phenomenon, 
but  this  could  not  be  substantiated  because  mercury  was  driven  into  the 
thermistor  when  pressure  was  applied  during  the  sample  run  with  the 
neoprene  sample.  This  small  change  in  volume  due  to  the  mercury  being 
driven  into  the  thermistor  was  enough  to  account  for  the  "volume  shift" 
observed  in  the  neoprene  sample  data.  It  was  known  that  mercury  was 
being  driven  into  the  thermistor  because  the  thermistor  inside  the 
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Fig.  10  —  Pycnometer  calibration  run  at  10° C.  Volume  (cm3),  shown 
on  the  ordinate,  represents  the  volume  of  the  mercury  in  the  bore  of 
the  acoustic  pycnometer.  A  sixth-degree  polynomial  is  shown  passing 
through  the  experimental  points. 
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Fig.  11  —  Pycnometer  calibration  run  at  25° C.  Volume  (cm3),  shown 
on  the  ordinate,  represents  the  volume  of  the  mercury  in  the  bore  of 
the  acoustic  pycnometer.  A  sixth-degree  polynomial  is  shown  passing 
through  the  experimental  points. 
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sample  chamber  short  circuited  during  pressurization  (i.e.,  exhibited 
zero  resistance),  and  when  the  pressure  was  released,  gave  a  temperature 
reading  similar  to  the  second  thermistor  on  the  outside  of  the  pycno¬ 
meter.  Moreover,  when  the  apparatus  was  disassembled,  small  droplets  of 
mercury  were  actually  found  to  be  embedded  in  the  epoxy  encasing  the 
thermistor.  This  behavior  indicated  that  there  was  a  flow  of  mercury 
into  and  out  of  the  thermistor  as  the  pressure  was  changed.  The 
thermistor  was  modified  so  mercury  could  not  short  It  out  during  the 
measurements  on  butyl-252. 

The  measured  values  of  V  made  in  the  sample  runs  are  given  in 

BUKh 

Figs.  12,  14,  and  17.  These  runs  show  no  evidence  of  hysteresis,  as 
stated  earlier.  In  fact,  it  made  no  difference  If  all  of  the  sample 
data  were  lumped  together  or  analyzed  cycle  by  cycle.  A  curve,  in  the 
form  of  a  sixth-degree  polynomial,  was  fitted  to  the  lumped  data  points. 

Figures  13,  15,  and  18  show  the  actual  sample  volumes  ^  as  a 

function  of  pressure  at  constant  temperature  for  butyl-252  and  type-W 
neoprene.  Figure  15  is  a  plot  of  the  volume  of  butyl-252  at  25°C  as  a 
function  of  pressure.  This  plot  is  of  particular  interest  because  it 
can  give  an  Idea  of  the  accuracy  of  the  volume  determination  when  using 
the  acoustic  pycnometer.  Theoretically,  this  zero-pressure  value  of  the 
volume  should  coincide  with  the  value  of  the  volume  of  the  butyl-252 
sample  determined  by  Archimedes'  principle.  The  two  determinations  do, 
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in  fact,  agree  closely.  The  difference  between  the  two  zero-pressure 
volume  values  of  the  butyl-252  is  on  the  order  of  0.22  percent.  Since 
the  zero-pressure  extrapolation  of  the  curve  fit  and  the  volume  deter¬ 
mination  by  Archimedes'  principle  are  in  such  close  agreement,  discarding 
data  below  2  MPa  appears  to  be  justified.  The  analysis  performed  to 
determine  V  for  the  butyl-252  sphere  was  also  applied  to  the 

measurements  made  on  the  type-W  neoprene  sample  at  25°C.  The  results 
are  shown  in  Fig.  18.  The  difference  between  the  zero-pressure  extra¬ 
polation  of  the  curve  in  Fig.  18  and  the  volume  measurement  of  the 
type-W  neoprene  sample  made  by  Archimedes'  principle  yields  an  error  of 
only  0.30  percent. 

Figure  13  illustrates  the  behavior  of  the  volume  of  the  butyl-252 
sample  at  10°C  as  a  function  of  pressure.  No  direct  comparison  of  the 
zero-pressure  extrapolation  of  the  curve  fit  in  Fig.  13  to  a  measurement 
of  volume  by  Archimedes'  principle  can  be  made.  In  determining  the 
volume  by  Archimedes'  principle,  one  requires  the  water  to  be 
approximately  the  same  temperature  as  the  surrounding  air.  If  a  temper¬ 
ature  gradient  occurs  in  the  water,  convection  currents  in  the  water 
will  be  initiated  and  a  false  weighing  of  the  sample  in  the  water  will 
result.  Therefore,  unless  the  water  temperature  can  be  kept  constant, 
Archimedes'  principle  will  produce  erroneous  results.  Unfortunately,  no 
apparatus  for  determining  volumes  by  weighing  the  samples  at  different 
temperatures  was  available.  An  estimate  can  be  made  to  compare  the 
zero-pressure  volume  extrapolation  of  the  butyl-252  sample  at  10°C  in 
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Fig.  12  —  Pycnometer  sample  run  of  butyl-252  at  10°  C.  Volume 
(cm3),  shown  on  the  ordinate,  is  the  volume  of  the  mercury  in  the  bore 
of  the  acoustic  pycnometer  with  a  sample  of  butyl-252  present  in  the 
sample  chamber.  A  sixth-degree  polynomial  is  shown  passing  through 
the  raw  data. 
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Fig.  13  —  Volume  of  the  butyl-252  sample  at  10°  C.  Actual  sample 
volume  is  given  as  a  function  of  pressure.  A  sixth-degree  polynomial 
is  shown  passing  through  the  points  calculated  from  Eq.  (lib). 
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Fig.  14  —  Pycnometer  sample  run  of  butyl-252  at  25°  C.  Volume 
(cm3 ),  shown  on  the  ordinate  is  the  volume  of  the  mercury  in  the  bore 
of  the  acoustic  pycnometer  with  a  sample  of  butyl-252  present  in  the 
sample  chamber.  A  sixth-degree  polynomial  is  shown  passing  through 
the  raw  data. 
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Fig.  15  —  Volume  of  the  butyl-252  sample  of  25°  C.  Actual  sample 
volume  is  given  as  a  function  of  pressure.  A  sixth-degree  polynomial 
is  shown  passing  through  the  points  calculated  from  Eq.  (lib). 
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Fig.  17  —  Pycnometer  sample  run  of  type-W  neoprene  at  25°  C. 
Volume  (cm3),  shown  on  the  ordinate,  is  the  volume  of  the  mercury 
in  the  bore  of  the  acoustic  pycnometer  with  a  sample  of  type-W  neo¬ 
prene  present  in  the  sample  chamber.  A  sixth-degree  polynomial  is 
shown  passing  through  the  raw  data. 
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Fig.  18  —  Volume  of  the  type-W  neoprene  sample  at  25°  C.  Actual 
sample  volume  is  given  as  a  function  of  pressure.  A  sixth-degree 
polynomial  is  shown. 
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Fig.  13  with  an  independently  determined  value  of  the  volume  of  the 
sample  by  using  the  value  of  the  butyl-252  sample  at  25°C  measured  by 
Archimedes'  principle  and  the  volume  coefficient  of  expansion  for  butyl 
rubber  [9],  The  value  used  for  the  coefficient  of  thermal  expansion  is 
57  x  10  ^/°C.  Performing  this  analysis,  one  calculates  a  value 
of  0.1  percent  for  the  difference  between  the  zero-pressure  extrapolation 
of  the  curve  fit  in  Fig.  13  and  the  volume  calculated  using  the  volume 
coefficient  of  thermal  expansion  for  butyl  rubber. 

Figure  16  shows  a  comparison  between  the  volume  of  the  butyl-252 
sample  at  10°C  and  25°C.  At  low  pressures  (less  than  5  MPa),  the  plot 
of  the  measurement  made  on  butyl-252  at  10°C  does  not  exhibit  the 
same  amount  of  curvature.  This  might  be  attributed  to  a  different 
compliance  of  the  rubber  at  the  lower  temperature. 

To  obtain  an  idea  of  the  actual  precision  of  the  acoustic  pycnometer 
measurement  technique,  the  observed  deviation  is  compared  to  the  random 

i 

deviation  computed  by  the  error  analysis  in  Appendix  B.  The  precision 
of  the  experiment  is  quite  good,  as  is  seen  ..n  the  graphs.  The  least- 
squares  curves  pass  directly  through  the  clusters  of  data  points.  The 
worst  case  of  deviation  is  observed  in  Fig.  15,  the  plot  for  the  sample 
volume  of  butyl-252  at  25°C.  However,  even  at  2  MPa,  the  worst  case, 
the  observed  deviation  is  only  0.17  percent.  The  maximum  observed 
deviation  for  the  plot  of  the  sample  volume  of  butyl-252  at  10°C  in 
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Fig.  12  is  0.09  percent.  The  calculated  random  deviation  when  the 
error  analysis  was  performed  is  on  the  order  of  0.1  percent. 

The  results  obtained  for  the  type-W  neoprene  sample  are  not 
presented  with  the  same  degree  of  confidence  as  the  results  for  the 
butyl-252  sample.  As  has  been  mentioned,  mercury  apparently  flowed 
into  and  out  of  the  thermistor  located  inside  the  sample  chamber  of  the 
acoustic  pycnometer  when  pressure  was  applied  and  released.  Since 
the  temperature  could  not  be  monitored  as  accurately  because  only  the 
thermistor  on  the  outside  of  the  pycnometer  was  functioning,  thermal 
equilibrium  was  not  assured.  Mercury  intrusion  into  the  thermistor 
in  the  sample  chamber  also  meant  that  a  small  amount  of  mercury  was 
being  lost  from  the  pycnometer  into  the  thermistor  inside  the  sample 
chamber,  so  the  volume  change  of  mercury  in  the  bore  was  not  due  just 
to  the  volume  change  in  the  sample.  This  thermistor  problem  is  also 
the  reason  why  no  data  were  obtained  on  the  type-W  neoprene  sample  at 
10°C  and  occasioned  the  design  of  the  thermistor  plug. 

The  sixth-degree-polynomial  least-squares  coefficients  for  the 
volumes  of  the  samples  of  butyl-252  and  type-W  neoprene  at  25°C  are 
respectively  given  in  Tables  I  and  11.  The  sixth-degree-polynomial 
coefficients  for  the  volume  of  the  butyl-252  sample  at  10 °c 
are  also  presented  in  Table  I. 

From  the  experimental  results,  the  isothermal  bulk  modulus  (see 
Appendix  H)  can  be  calculated  by  Eq.  (H3).  The  results  of  the  bulk- 
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Table  I  —  Least-squares  coefficients  used  for  the  calculation  of  the 
volume  of  the  butyl-252  sample  at  10°  C  and  25°  C.  The  equation  for 
the  volume  of  the  sample  has  the  form:  VS(P)  =  Aj  +  A2P  +  A3P2 
+  A4P3  +  A5P4  +  A6P5  +  A7P6. 


Butyl 

-252  Least-Squares  Coefficients,  P 

=  Ambient  to  70  MPa 

T  =  10°C 

T  =  25°C 

A1 

+ 

0.881761  x  101 

+ 

CU886270  x  101 

A2 

- 

0.439446  x  10-3 

- 

0.106370  x  10“  2 

A3 

- 

0.175415  x  10-3 

- 

0.157199  x  10“  3 

\ 

+ 

0.718107  x  10-5 

+ 

0.633567  x  10“5 

A5 

- 

0.147986  x  10-6 

- 

0.119845  x  10“6 

-  A6 

+ 

0.152569  x  10“8 

+ 

0.108048  x  10“8 

A7 

- 

0.625368  x  10“ 11 

0.367876  x  10"11 

Table  II  —  Least-squares  coefficients  used  for  the 
calculation  of  the  volume  of  the  type-W  neoprene 
sample  at  25°  C. 


Type-W  Neoprene  Coefficients 

P  =  Ambient  to  70  MPa,  T  =  25°C 

A1 

+ 

0.8803583  x  101 

A2 

- 

0.1917722  x  10"1 

A3 

+ 

0.8394465  x  10-3 

A4 

- 

0.1701149  x  10“4 

A5 

+ 

0.7423981  x  10"7 

A6 

+ 

0.1733432  x  10“8 

a7 

- 

0.1609281  x  10" 10 

modulus  calculation  are  given  in  Fig.  Hi.  As  shown  in  Fig.  HI,  the 
isothermal  bulk  modulus  of  the  butyl-252  sample  rapidly  decreases  upon 
initial  compression  and  then  becomes  relatively  constant  at  pressures 
greater  than  10  MPa.  The  bulk-moduli  curves  of  the  butyl-252  sample 
intersect  at  high  pressure  (~65  MPa).  This  is  probably  a  numerical 
artifact  associated  with  taking  the  derivatives  of  the  sample-volume 
curves  in  Fig.  16,  which  are  needed  for  the  calculation  of  the  bulk 
modulus.  The  drastic  change  in  the  low-pressure  (less  than  10  MPa) 
bulk-moduli  results  In  Fi.g.  'll  might  be  due  to  some  sort  of  rearrangement 
of  the  molecular  structure  of  the  butyl-252  elastomer  sample  on  a  micro¬ 
scopic  scale.  At  this  time,  this  conclusion  is  purely  speculative.  At 
pressures  greater  than  10  MPa,  the  bulk-modulus  of  butyl-252  increases 
slightly  as  a  function  of  pressure,  or  the  rubber  becomes  stiffer  with 
increasing  pressure  [10].  This  behavior  is  also  shown  in  Fig.  HI. 

VI.  SUMMARY  AND  CONCLUSIONS 

The  experimental  results  given  in  the  preceding  section  demonstrate 
that  the  acoustic  pycnometer  is  an  instrument  for  precise  determination 
of  the  volume  of  an  elastomer  sample  as  a  function  of  pressure  at  any 
temperature. 

The  only  measure  presently  available  of  the  accuracy  of  the  pycno¬ 
meter  results  is  obtained  by  comparing  the  zero-pressure  extrapolation 
of  the  least-squares  curve  fit  of  the  volume  of  the  samples  to  the 
volume  of  the  samples  when  measured  by  Archimedes’  principle  at  25°C. 
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Since  this  comparison  yielded  errors  of  0.22  percent  and  0.30  percent, 
respectively,  for  the  samples  of  butyl-252  and  type-W  neoprene,  it  can 
be  assumed  the  method  used  in  this  research  is  accurate  as  well  as 
precise. 
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APPENDIX  A 


PULSE-ECHO  OVERLAP  TECHNIQUE 

The  time-of -flight  measurement  described  earlier  is  done  by  a 
method  called  the  pulse-echo  overlap  technique  [I].  This  technique  was 
originally  designed  to  accurately  measure  the  speed  of  sound  in  a 
material.  One  calculates  the  speed  of  sound  in  a  material  from  measure¬ 
ments  of  the  time  it  takes  for  an  ultrasonic  pulse  to  propagate  through 
the  material  and  from  knowledge  of  the  thickness  of  the  material  sample. 

Instead  of  using  the  pulse-echo  overlap  technique  to  measure  the 
unknown  speed  cf  sound  in  a  sample  of  material  by  measuring  the  time  of 
flight  in  the  sample,  it  was  used  in  this  investigation  to  determine  the 
height  of  a  column  of  mercury  (the  sound  speed  in  which  is  known)  in  the 
bore  of  the  pycnometer  by  Eq.  (8).  In  the  pulse-echo  overlap  method  (Fig. 
Al)  the  cw  oscillator  frequency  is  adjusted  so  that  its  period  is  exactly 
equal  to  the  round-trip  time  in  the  sample  [12].  A  standard  oscilloscope 
is  used  to  observe  the  pulse-echo  train  shown  in  Fig.  5.  The  time 
difference  between  two  consecutive  echoes  is  the  time  it  takes  for  a 
sound  pulse  to  travel  from  the  transducer  in  the  base  of  the  bore  of  the 
pycnometer  through  the  mercury  column,  to  reflect  from  the  float,  and  to 
return  to  the  transducer.  The  Matec  122-B  dual  decade-delay  generator, 
shown  in  Fig.  Al ,  uses  strobe  signals  (displayed  via  the  oscilloscope) 
to  highlight  the  first  and  second  echoes  in  the  pulse-echo  train.  The 
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Fig.  A1  —  Schematic  of  acoustic-pycnometer  measurement  system 
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Hewlett-Packard  frequency  synthesizer  Is  set  approximately  to  the  time 
difference  between  the  first  and  second  echoes  as  determined  by  a  rough 


time  measurement  using  the  oscilloscope.  These  two  echoes  can  be  over¬ 
lapped  or  superimposed  when  the  horizontal  sweep  of  the  oscilloscope  Is 
triggered  by  the  cw  triggers  In  the  Model  122-B  (Sweep  Sync  Out-Direct). 
By  carefully  adjusting  the  vernier  on  the  frequency  synthesizer  and  by 
adjusting  the  time  base  on  the  oscilloscope,  a  precise  overlap  can  be 
achieved  (see  Fig.  A2);  hence  an  accurate  measurement  is  made  of  the 
time-of-f light  in  terms  of  a  precisely  set  synthesizer  frequency.  When 
aligning  the  echoes  on  the  oscilloscope,  the  earliest  distinguishable 
feature  is  used  as  a  reference  to  avoid  mismatching  the  echoes.  The 
only  constraint  is  that  the  total  sweep  time  must  be  less  than  the 
round-trip  time  in  the  mercury  column  to  allow  for  retrace  time  in  the 
oscilloscope  [13]. 
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Fig.  A2  —  Pulse-echo  overlap:  timebase  =  0.5  jus.  The 
earliest  arriving  echo  has  the  greatest  amplitude. 


APPENDIX  B 


ERROR  ANALYSIS  OF  THE  ACOUSTIC  PYCNOMETER  MEASUREMENT 


The  precision  of  the  method  used  in  the  acoustic-pycnometer 
measurement  can  be  investigated  by  examining  the  random  deviation  of  the 
components  of  Eq.  (14b).  The  maximum  resultant  effect  due  to  several 
deviations  in  any  measurement  is  simply  the  arithmetic  sum  of  the 
individual  deviations  in  the  measurement  [14],  or 


RES 


+  A 


y 


(Bl) 


where 


Ares  *  deviation  the  resultant  effect 
through  Aj=  the  deviations  attributed  to  the  various  components 
of  the  measurement. 

Each  deviation  in  Eq.  (Bl)  can  be  expressed  as 


j 


3M  . 
9*j  3 


(B2> 


where 

M  *  the  formula  used  to  calculate  the  quantity  desired  and  is 
a  function  of  the  properties  measured 
nij  *  the  property  in  which  M  exhibits  a  deviation 
6^  *  the  actual  deviation  in  a  particular  property  measured. 
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Equation  (Bl)  represents  a  maximum  possible  random  error  because  the 
sums  of  the  individual  deviations  are  taken  to  be  positive  quantities 
but  in  reality  could  have  different  signs.  The  most  probable  resultant 
deviation  A^g,  according  to  method  of  least-squares  [15],  is 


^  +  A^  +  A  ^  + 

L  2  3  *  * 


.  .  +  A 


An  expression  for  the  resultant  deviation  of  the  calculation  of  the 
volume  of  the  sample  in  the  pycnometer  measurement  is  obtained  by  applying 
Eq.  (Bl)  to  Eq.  (14b).  The  resultant  deviation  is 


4VS»MPLE(P'T>  -  “sAMPLe'W  +  ®BORE  +  “BORE  +  4V(P-T)-  (B4) 


The  individual  deviations,  in  Eq.  (B4),  are  in  turn  the  sum  of 

component  deviations  experienced  when  making  their  measurement. 

Therefore,  the  individual  deviation  of  V„.,_-_(P  ,T  )  is 

SAHrLfc  O  O 

4Vsrmple(po’V  ■  {4m.h  (4VoJ +  {4oh,o}  ■ 


(B5a) 


where 


AVSAMPLE(Po,TV  “  VSAMPLE(ma'  "l^O’  pH20(Po’To^* 


The  expression 


ha} 


refers  to  the  deviation  in  volume  due  to 


weighing  the  mass  of  the  elastomer  sample  in  air.  The  notation  used 


to  express  the  deviation  of  the  mass  of  the  elastomer  sample  in  air 

W is 


8V 

SAMPLE 

3m 

a 


(B5b) 


The  individual  deviation  associated  with  measuring  the  volume  of 
the  mercury  in  the  bore  is 

4vbore  ’{4rb}+  H  +  <B6*> 


since 


VBORE  "  VBORE  (rb(p»T>»  CHg(P,T)*  f)*  (B6b) 

The  expression  for  the  individual  deviation  for  the  last  term  in  Eq. 

(B4)  is 

SV(P.T)  -  {am'}  *j«.j  ♦  {4pHg}  +  {4PSAMPBE^Po'Po^} 

+  {4Po}+  {4To}  +  H  +{4T1-  <B7> 

where 

V<P-T>  -»<•’.  ».  »„,<*„>•  VSAMpLE(P0I0),  P0,T0,P,T). 

The  individual  deviations  given  in  Eqs.  (B5a),  (B6a),  and  (B7) 


can  be  analyzed  by  using  Eq.  (B2).  Applying  Eq.  (B2)  to  Eq.  (B5a),  the 
deviation  of  the  initial  sample  volume,  one  obtains 


AVSA>!PLE^Po,To) 


8V 

SAMPLE  , 

4_ 

Dv 

SAMPLE  r 

am 

a 

T 

3%0  \° 

SAMPLE 

3pH20 


(B8a) 


After  taking  the  appropriate  derivatives,  one  finds  that  the 
equation  for  the  random  deviation  for  the  initial  sample  volume  is 


4We< VV  -  *■. +  ^  *V 

("a  +  V  a°H,0 

+  —2 - TT-  «. 


(B8b) 


In  Eq.  (B8b),  6ma,  Q,  and  6T  are  observed  deviations  associated  with 

measuring  the  mass  of  the  sample  in  air,  the  mass  of  the  sample  in 

water,  and  the  temperature,  respectively.  It  should  be  noted  that  an 

appropriate  estimate  Is  needed  for  Ap  ../AT.  This  estimate  is  obtained 

H2° 

by  using  the  table  of  the  density  of  water  p  _  as  a  function  of  temper- 
ature  and  taking  a  finite  difference  for  Ap„  n  between  a  temperature 
above  and  below  the  temperature  T  . 


By  the  preceding  analysis,  the  random  deviations  can  be  calculated 

I 

for  the  case  of  the  volumes  of  mercury  in  the  bore  AV  and  AV_„_  and 

BORE  BORE 


the  volume  of  mercury  associated  with  filling  the  pycnometer  to  a 


different  level  V(P,T).  Substituting  this  value  obtained  into  Eq. 
(B4),  one  can  obtain  an  estimate  of  the  resultant  deviation  for  the 
measurement  of  volume  of  any  sample  as  function  of  pressure  and 
temperature  in  the  acoustic  pycnometer. 


It  is  found  that  the  deviation  calculated  from  Eq.  (B4)  is 

3 

0.0088  cm  .  This  represents  a  random  deviation  in  V„ „(P,T)  of  0.1 

SAMrLn, 

percent.  This  estimate  of  0.1  percent  is  calculated  from  the  following 


data: 


P  =  0.0  MPa 


6P  =  ±0.1  MPa 


P  =  0.0  MPa 
o 


SP  =  ±0.1  MPa 
o 


T  =  25°C 


6T  =  ±0 . 2  °C 


T  =  25°C 
o 


6T  =  ±0.2°C 
o 


f  =  13.0  kHz 


6f  -  ±5.0  Hz 


6m  =  ±0.0003  g 
a 

Sm^,  Q  =  ±0.0003  g 
r,  =  0.5545  cm 

D 


cu  =  147150  cm/s 


Oy  »  3.6  X  10  /°C 


6^  *  0.0/MPa. 


APPENDIX  C 


f 


CALCULATION  OF  THE  VOLUME  OF  MERCURY  IN  THE 
BORE  OF  THE  ACOUSTIC  PYCNOMETER 


The  volume  of  the  mercury  in  the  bore  of  the  acoustic  pycnometer 
is  given  by 


V 


BORE 


AL  , 


(Cl) 


where 

A  =  area  of  the  bore 

L  =  height  of  the  mercury  in  the  bore. 


The  area  of  the  bore  in  Eq.  (Cl)  is  a  function  of  temperature  and 
pressure  governed  by  the  coefficient  of  thermal  expansion  and  the 
coefficient  of  compressibility  of  the  pycnometer  container.  The  volume 
coefficient  of  thermal  expansion  is  given  by 


a 


V 


(C2) 


Therefore,  at  constant  pressure  the  change  In  volume  per  unit  volume  is. 


dV 

V 


UydT  . 


(C3) 


The  isothermal  volume  coefficient  of  compressibility  is 


Therefore,  at  constant  temperature  the  change  in  volume  per  unit  volume 


is 


dV 

V 


-  e.,dp 


(C5) 


The  expression  for  the  change  in  area  per  unit  area  as  a  function  of 
temperature  and  pressure  would  then  take  on  the  form 


dA 

A 


3  |_  \ 


dT 


-  B  dP 
V 


(C6) 


Integrating  Eq.  (C6),  one  obtains 

A(P,T)  =  irr^  exp[2(avT  -  ByP  )/3] . 


(C7) 


The  volume  VB0RE  of  the  mercury  in  the  bore  can  be  calculated  by 
Eq.  (Cl)  using  Eq.  (8)  and  Eq.  (C7).  It  is  found  to  be 


WP>T)  "  2  tfirrb  exP f “ ( ayT  -  BvP)/3] 


(C8) 


A  Hewlett-Packard  Model  9825  Calculator  is  used  to  calculate  V 

BORE 

from  Eq*  (C8)*  The  program  is  shown  on  the  following  page. 


(Sh 


0: 


8tart":fxd  5 


1:  enp  "Temp. 
[Celcius] " ,T 

2:  enp  "Pressure 
[MPa] " ,P 

3:  enp  "Freq. 
[kHz ] ",F 

4:  1460-*r6;  ,46->r 

7;  .  21-*r8 

5:  CKr9;  9e-6-+rlO 

6:  .  5545-+R 

7:  (r6-r7*T+r8* 

P)*100->C 

8:  prt  "CHg  [cm/ 
s]",C 

9:  .66667  (r9*T- 

rlO*P)+D 

10:  exp(D)-+E 

11:  .5*ttR+2*E*C/ 
(F*1000>+B 

12:  prt  "VHBG 
[ml]  »",B;spc  2 

13:  gto  2 

*1838 
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APPENDIX  D 


LEAST-SQUARES  ANALYSIS  AND  COMPUTER  PROGRAM 

In  order  to  use  Eq.  (14b)  to  determine  the  volume  of  an  elastomer 
sample  as  a  function  temperature  and  pressure,  it  is  necessary  to  de¬ 
termine  an  expression  for  the  volume  of  the  mercury  in  the  bore  of  the 
pycnometer  during  the  sample  run  and  the  calibration  run.  This  can  be 
accomplished  by  obtaining  an  equation  of  "best  fit"  to  the  experimental 
data  by  the  method  of  least  squares. 


Since  the  experimental  data  did  not  follow  a  linear  relation 
when  plotted,  a  general  least-squares  solution  is  employed  in  the  form 
of  a  rational  fraction: 


y(x) 


2  3  A  5  6 

A,  +  A_x  +  A,x  +  A. x  +  Acx  +  A,x  +  A,x 
1  2 _ 3 _ 4 _ 5 _ 6 _ 7 _ 

2  3  A  5  6  * 

A8+A9X+A10X  +A11X  +A12x  +A13x  +  A14x 


(Dl) 


In  Eq.  (Dl),  the  quantities  x  and  y  are  the  experimental  observations. 
The  coefficients  A^  through  A^  are  determined  by  the  method  of  least 
squares  as  a  best  fit  to  the  observed  values  of  x  and  y.  Rearranging 
Eq.  (Dl)  into  a  more  convenient  form,  one  obtains 


AQy  +  AQxy  +  Air)x  y  +  A,,x  y  +  A,0x  v  +  Anx  y  +  A,ax  v 


*11 


12' 


'13' 


'14' 


2  3  4  3b 

A.  -  A?x  ~  A-,x  -  A^x  -  A t- x  -  A^x~  -  A^x  =  0. 


(D2) 


Equation  (D2)  represents  an  equation  of  best  fit  for  a  sixth-order 
rational  fraction  if  no  experimental  errors  are  present  in  the 
observed  quantities  x  and  y.  However,  this  is  never  the  case,  since 
all  measurements  are  subject  to  error.  Suppose  one  has  n  pairs  of  ex¬ 
perimental  measurements  x^,  y  (i  *  1,  2,  .  .  .  n).  Then  Eq.  (D2) 
can  be  written  in  terms  of  the  observed  quantities  x^  and  y^  and  the 
experimental  error  associated  with  measuring 

Agyi  +  AgX^  +  AjqX^j  + 

A8y2  +  A9X2y2  +  AlOx2y2  +••••- 

Vi  +  Viyi  +  AioVi  +  •  •  •  •  - 


x^  and  y^ 

.  6 

A7x1  “  V1 

.  6  _ 

A^2  v2 

A7xi  =  vi  (D3) 


A8yn  +  A9Vn  +  A 


vlO*n  yn 


-  A7x^ 
7  n 


In  Eq.  (D3),  the  subscripts  associated  with  the  quantities  x^,  y^, 
and  refer  to  the  i-th  observation. 


The  least-squares  principle  states  that  the  coefficients  A1 
through  A  in  Eq.  (D3),  are  the  best  fit  to  the  data  when  the  sum 
of  the  squares  of  the  errors  is  a  minimum  [16],  or 
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0 


2> 


i  =  i 


3v 


n 


0 


(D5) 


9v 


9A 


14 


=  0. 


Equation  (D5)  produces  a  system  of  linear  simultaneous  equations  in 
terms  of  the  observed  quantities  x^,  y^.  In  order  that  this  system 
of  simultaneous  equations  is  not  linearly  dependent,  one  of  the 
coefficients  must  be  set  to  a  known  value.  To  solve  a  system  of 
simultaneous  equations  of  this  nature,  in  order  to  obtain  the  coeffi¬ 
cients  Aj  through  A^,  a  digital  computer  is  needed.  Equation  (D5) 
can  be  written  in  matrix  form  as 

[E]  [A]  -  0.  (D6) 

In  Eq.  (D6),  the  large  matrix  [E]  represents  the  terms  in  Eq. 

(D5)  after  the  appropriate  derivatives  have  been  taken  and  multiplied 


fc. 

iAi 
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by  the  error  term  v.  The  [A]  matrix  is  a  column  matrix  that  represents 
the  coefficients  which  best  fit  the  observed  ata. 


The  computer  program,  RAGT.FTN,  using  NRL-USRD's  PDP-11/45 
computer,  was  written  in  FORTRAN  to  do  the  least-squares  analysis.  The 
program  was  written  in  a  general  fashion  to  try  and  accommodate  various 
types  of  data.  For  instance,  a  polynomial  curve  fit  is  obtained  by 
simply  initializing  the  coefficient  Ag  to  1  and  by  setting  coefficients 
kg  through  A^  to  zero.  If  a  lower-degree  polynomial  fit  is  desired, 
one  simply  suppresses  the  appropriate  coefficient  by  setting  it  to 
zero.  The  subroutine  SELECT. FTN  allows  the  user  to  select  any 
coefficient  for  initialization.  The  subroutine  SUPRES.FTN  sets  any 
desired  coefficient  to  zero  to  lower  the  degree  of  polynomial.  The 
subroutine  DSIMEQ.FTN  solves  the  system  of  simultaneous  equations 
shown  in  Eq.  (D6),  by  means  of  Gaussian  elimination  in  double  precision 
Once  the  least-squares  coefficients,  represented  by  the  column  matrix 
[A]  in  Eq.  (D6),  are  calculated,  the  subroutine  YCALCU.FTN,  uses  these 
least-squares  coefficients  to  recalculate  the  function  y(x)  in  Eq.  (Dl) 
A  listing  of  RAGT.FTN  and  all  the  subroutines  used  follows: 
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C  PROGRAM  DIRECTORY 

C  DATA  TO  DE  ENTERED 

C  M=*  OF  DATA  POINTS 

C  XPRESS! 100)=  ARRAY  TO  HOLD  PRESSURE  DATA 

C  YVOL< 100 )  =  ARRAY  TO  HOLD  VHGB  DATA 

C  A! 14. 14)=  ARRAY  TO  HOLD  VARIOUS  SUMS 


C  B!N)  =  COLUMN  MATRIX 

C  N=  NUMBER  OF  SIMULTANEOUS  EQUATIONS 

BYTE  ANSC4 ) 

DIMENSION  XPRESSC 100) ,YVOL( 100) ,XINC( 100) .YCALC(IOO) 

DOUBLE  PRECISION  A! 14 . 14) ,B! 14) ,X!14> , VALUE! 14) ,DX1 .DX2.DX3, 
1 DX4 , DX5 , DX6 , DX7 , DXB , DX9 , DX1 0 , DX1 1 .DXI2.DY1 .DY2.DXIY1 ,DX1Y2, 
IDX2Y1 , DX2V2 , DX3Y1 .DX3Y2.DX4Y1 ,DX4Y2,DX5Y1 , DX5Y2.DX6YI , DX6Y2 
IDX7Y1 , DX7Y2 , DXBY 1 , DXHY2 , DX9 Y 1 , DX9Y2 , DX10YI .DXI0Y2.DX1 1Y1 , 
1DX1 IY2.DX12Y1 .DX12Y2.X1 ,X2 . X3 ,X4 ,X5 , X6 ,X7 ,X8,X9 ,X10 ,X1 1 .X12 
1Y1 .Y2.X1Y1 .X1Y2.X2Y1 .X2Y2.X3Y1 .X3Y2.X4Y1 .X4Y2.X5Y1 ,X5Y2, 
1X6Y1 .X6Y2.X7Y1 .X7Y2.X0Y1 .XHY2.X9Y1 .X9Y2.X10Y1 .X10Y2.X11Y1 , 
1X1 1Y2.X12Y1 , XI 2Y2.DUM, SCALE 
WRITE! 3. 1000) 

1000  FORMAT! x,’ SUSE  OLD  DATA?*) 

READ! 5 , 1010 , END=800 ) ANS 
1010  FORMAT (4A1 ) 

I F! ANS! 1 ) . Eft. ' N* )C0  TO  870 

OPEN! UN1T= 1 , NAME= ’  RAGT.DAT’ ,TYPE= ’ OLD* ) 

READ! 1 . 1Q)M 
DO  860  1=1 ,M 

READ! 1 , 12 )XPRESS ! I ) .YVOL! I ) 

860  CONTINUE 

GO  TO  110 

870  OPEN ! UN IT=  1  , NAIIE=  ’  RACT. DAT’  ,TYPE=  ’  NEW  ) 

C  ALLOCATE  STORAGE 

DO  10  1=1,100 
XPRESS! I ) =0 .0 
YVOL! I ) =0 .0 
10  CONTINUE 

WRITE! 5,15) 

15  FORMAT! 5X, ’  $ ENTER  THE  NUMBER  OF  DATA  POINTS*) 

READ! 3 , 18)M 

10  FORMAT! 1 1 0 ) 

DO  100  .1=1,  N 
WRITE! 5, 1 1 ) 

1 1  FORMAT! IX, '  XPRESS , YVOL’ ) 

READ! 5, 121XPRESS! J> , YVOL! J) 

12  FORM AT !2F 10.0) 

100  CONTINUE 

Will  IE!  I  .  IU)M 
DO  40  1=1 ,M 

WRITE! I , 1020) XPRESS! I ) ,YV0L!1 ) 

1020  FORMAT! F 1 0 . 5 , F 1 0 . 5 ) 

•10  CONTINUE 

110  CLOSE! UNIT= 1 ) 

WRITE! 5. 1040) 

1040  FORMAT!  /* , '  $ ENTER  HIE  NUMBER  OF  SIMULTANEOUS  EQUATIONS’ ) 
READ! 5, 1 050 ) N 
lOUO  FORMAT! 14) 

DATA  XI  .X2.X3.X4,Xr».X6.X7,X8,X9,X!0,Xl  I  .X12.Y1  .Y2.X1Y1  ,X1Y2. 
1X2Y1 ,X2Y2,X3YI .X3Y2.X4Y1 ,X4Y2,X5Y1 ,X3Y2,X6Y1 .X6Y2.X7YI . 
1X7Y2.XHYI .X8Y2.X9Y1 .X9Y2.X10Y1 .X10Y2.XI 1Y1 ,X1 1Y2.XI2Y1 , 

I  XI  2Y2/0U  M) .  0U+02/ 

DO  50  1=1 ,M 

C  SET  UP  DUMMY  STORAGE 

DXI=  DIJI.E!  Xl’RESS!  I  )  ) 

RX2  =  I)XI*UX1 
DX3=DXICDX2 


0X4 =0X2* 0X2 
DX5=DX2*DX3 
0X6=DX3*0X3 
0X7  =  0X1*11X6 
0X0=0X2*11X6 
0X9  =  0X3*11X6 
0X10=0X4*0X6 
I)X1  1  =0X5*UX6 
0X12=0X6*0X6 

dyi=obu:(Yvol(i>> 

0Y2=0Y1*DY1 

I1X1Y1=0X1*DY1 

0X1 Y2=0X1 *0Y2 

1)X2Y  1  =DX2*DY  1 

0X2  Y2  =  11X2*0  Y2 

l)X3Yl  =DX3*DY1 

0X3Y2=0X3*0Y2 

UX4Y1  =DX4*0Y1 

DX4Y2=OX4*OY2 

0X5Y1=0X5*0Y1 

0X5Y2=DX5*0Y2 

DX6Y1 =DX6*0Y  1 

DX6Y2=DX6*DY2 

DX7Y1=DX7*DY1 

DX7Y2=DX7*0Y2 

DXHY 1 =  DXB*OY 1 

DX0Y2=DXn*DY2 

DX9Y1 =DX9*DYl 

0X9Y2=DX9*0Y2 

DXI0Y1 =oxi o*oyi 

DX10Y2  =  DX10*I1Y2 
DXl 1Y1 =0X1 1 *DY1 
OX I 1Y2=0X1 !*0Y2 
0X1 2Y1 =DX12*0Y1 
0X12Y2=DX12*DY2 

start  accumulating  SOWS 

XI =X1 +0X1 

X2=X2+0X2 

X3=X3+0X3 

X4=X4+0X4 

X3=X5+DX5 

X6=X6+0X6 

X7=X7+0X7 

xa=xa+oxo 

X9=X9+0X9 
X10=X10+DXI 0 

XI  1 =X1 1+0X1 l 
X12=X!2+0X12 
Y!=YI+0Y1 
Y2= Y2+DY2 
XI Y1 =X1 Y1 +0X1 Yl 
XI Y2=XI Y2+0X1Y2 
X2Y I =X2Y1 +0X2Y 1 
X2Y2= X2Y2+ 0X2Y2 
X3Y1 =X3Y1 +DX3Y1 
X3Y2=X3Y2+0X3Y2 
X4  Y I  =  X4  Y  l  +  0X4Y 1 
X4Y2=X4Y2+0X4Y2 
X5YI =X5Y1 +DX5YI 
X5Y2=X5Y2+0X5Y2 
X6Y1 =X6YI +0X6Y1 
X6Y2=X6Y2+0X6Y2 
X7Y1 =X7YI +0X7Y 1 
X7Y2  =  X7Y2+I*X7Y2 
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XI)  Yl  =XIIY1  +UX1JYI 
XHY2  =  XIIY2+I>XIIY2 
X‘>Y I  =  X9Y  1  +  I)X‘)Y  I 
X*>Y2  =  X‘)Y2fDX9Y2 
XI0Y1 =XIOYI +DXIOY1 
XIOY2=XI 0Y2  +  IIX1 0Y2 
XI  I Y l =  X I  IYI+DXI IY1 
XI  1 Y2=XI 1Y2  +  DXI IY2 
X12Y1  =X12YI +11XI2YI 
XI 2Y2  =  X I 2Y2  HIXI2Y2 
CONTINUE 

STONE  SUMS  IN  MAT1UCES  ACN.N)  AND  3(N) 

A  (  1  , 1 )  =  Y2 

A(2 , 1 >=X1 Y2 

A (2 ,2 ) =X2Y2 

A Cl.  1  >=X2Y2 

AC)  ,2)  =  X3Y2 

A  C) ,  3  )  =  X4Y2 

A<4, 1 )=X3Y2 

A ( 4 , 2) =X4Y2 

A ( 4 ,3 )  =  Xi5Y2 

A ( 4 , 4 ) -X6Y2 

A(  5 , I  )=X4Y2 

A  (5, 2) =X3Y2 

A  <  5 , 3 ) -X6Y2 

A ( 3 , 4 ) =  X7Y2 

A<r>.r»)=xav2 

A  <  6 , 1 ) =X5Y2 
A(6,2)=X6Y2 
A (0 , 3 ) -X7Y2 
AC  f>  ,4  >  =  X8Y2 
A(6.5)=X*>Y2 
A( 6 , 6 ) =XI 0Y2 
A (7, 1 )=X6Y2 
A( 7 , 2 ) =X7 Y2 
A  ( 7 , 3  )  =  XI)Y2 
A (7 , 4 ) -X9Y2 
A  <  7 . 5 )  =  X I 0Y2 
A  ( 7 , 6  )  =  X  I  I Y2 
A ( 7 . 7 ) =X I 2Y2 
A  (II.  1  )  =  -Y 1 
A  Cl . 2 )  =  -X I Yl 
A(1I,3>  =  -X2YI 
AC  II  ,4 )  =  -X3Y I 
A  (II .  5 )  =  -X4Y I 
A(II.6>  =  -X5Y1 
A  (II  ,7  )  -  -X6Y I 
ACI.Il)=n 
AC),  1  »  =  -XI  Yl 
A  < ')  .  2  )  -  -a'2Y  I 
AC),3)-~X3YI 
A('l,4)=  -X4YI 
A  0),._»)  =  -X5Yl 
A  ( •>  .  (. )  =  -Xu  Y  I 
A  l  ‘1 . 7  i  -  -\T't  I 
\  < •> .  •  l » -  X  I 
\<  •>  .1  •  -.2 
A  i  10  .  I  '  -  l'i  I 
\  1  I  ti  .  2  i  -  ' ■ : l Y  ! 

\  <  in  .3  I  -  -  .  I  ■■  ! 


A  (  1 0  .  IJ  >  =X2 
A(  io,li)-x:i 
A( 10. I0)=X4 
A(  1  1  ,  1  )=-X3Yl 
A(  I  I  ,  2  )  =  -  X4  Y 1 
A  < I 1  ,3>=-X3YI 
A(ll . 4 ) =  -X6Y I 
Adi  ,  5  )  =-X7YI 
A  (  d  ,  6  )  =  -XHY l 
A  (  1  1  ,7  )  =  -X9Y  1 
A ( 11 , 0 ) =  X3 
Adi  ,<>)=:i4 
A ( 11 , 1 0 )  =  X 3 
A ( I  I , 1 1 )  =  XA 
A (12, 1 )=-X4Yl 
A ( 12,2) =-X3Yl 
A  < 12,3 ) --X6Y1 
A  (  12,4)  = -X7Y 1 
A  (  12,3)  =  -XHY  1 
A  (  1 2 ,  A  )  =  -X9Y  1 
A( 12,7) =-Xl OYl 
A( 12 ,U)-X4 
A ( 12,9) =X3 
A( 12, 10)=X6 
A( 12. 1 1 )=X7 
A( 12, 12)=XH 
A( 13, 1 ) = -X3Y1 
A  <  13,2)  =  -X6Y1 
A  ( 1 3 , 3 )  =  -X7Y1 
A  < 13,4 )  =  -XllY 1 
A (  13.3) =  -X9Y 1 
A ( 13,6 )  =-Xl OYl 
A ( I  3,7)  =  -XI  1Y1 
A (  I3,H)=X3 
A ( 13,9 )  =  X6 
A ( 13, I0)=X7 
A  (  1  3 , 11  )  -Xll 
A( 13, 1 2 ) =X9 
A( 13, 13 ) =X1 0 
A (14, 1 )=-X6Yl 
A( 14,2) - -X7Y 1 
A ( 14,3)=-X»YI 
A ( 14,4) --X9Y 1 
A ( I  4 , 5  )  =  -X 1 OYl 
A ( 1 4 , 6  > --XI I Y 1 
A ( 14,7) =-Xl 2Y1 
A ( 14  ,H) =  X6 
A ( 14,9)= X7 

A ( 1 4 , 1 0 )  =  xa 

A ( 1 4 , 11 ) =  X9 
A  (  I  4 . 1  2  )  -  X  I  0 
A ( 1 4 ,  I  3 ) -X 11 
A (  14, I  4  )  -  X I  2 

HUE  TO  Till'.  SYNMKTKY  OK  MATRIX  A  ( N ,  W)  ,  THE  RKMAINDEU  OF  THE 
riAIUIX  IAN  STORED  MY  TIIE  USE  OF  A  l.OOF  US  I  NO 
1  III.  Mil  .AT  ION  A(K  ,.n  A(.I,K) 

MO  121)  ).  2  .  M 

:■  kki 

.  'II  ||  3  .1  1  .  ‘■’K 

\  .  .1  .  )  '  M  ¥  .  .1  ) 

'■.".i  l  l  HIM' 

■  ■•ni  l  M  ! 

i>K  I  II!'  SOM  I  HO  I  ACTOR 


1 1  ■ 


on 


1220 

1230 

C 


70 

75 


1200 


D 

D1060 

C 

1300 

C 

COO 

2000 

«ioo 


FORMAT !  / ,  ’  SS  ENTER  THE  SCALING  FACTOR’) 

READI5, 1230)SCAI.E 
FORMAT! D24 . 17) 

scale  'i he  matrix  acn.N) 

DUM  =  (0  .  lOOOD+Ol  ) /SCALE 
DO  75  J  =  1 ,N 
1)0  70  K=  1  ,  N 

A! J,K)=A!J,K)*DUM 

CONTINUE 

CONTINUE 

T11E  SUIIHOUTINE  SELECT  ALLOWS  YOU  TO  SELECT  ANY  OF  THE  CO¬ 
EFFICIENTS  FOR  INITIALIZATION 
CALL  SELECT (N, A. B, VALUE) 

WRITE! 5. I200) 

FORMAT!/,’  ^SUPPRESS  ANY  COEFFICIENT?’) 

READ! 5, 1010, END =800 )ANS 
IF!ANS! I ).EQ.’N’ )GO  TO  1300 
CALL  SUPRES! N , A , B,X) 

WRITEI5, 1060) 

FORMAT!’  READY  TO  CALL  DSIMEQ’ ) 

CALL  SUBROUTINE  TO  SOLVE  SIMULTAHEOUS  EQ. 

CALL  DSIMECKN.A.B.X) 

CALL  SUBROUTINE  YCAL  TO  CALCULATE  VALUES  OF  YVOL 
CALL  YCAL! N ,X , YCALC ,XINC) 

WRITE! 5 ,2000 ) 

FORMAT! ’  PROGRAM  EXITS! ’ ) 

CO  TO  900 
CALL  EXIT 
END 
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SUBROUTINE  DSINEOIN.A.B.X) 

DOUBLE  PRECIS  1  ON  A(  1 4 , 14 )  , B<  1 4 )  ,X(1 4 )  , DET. SUN ,C 
DIMENSION  Jl’RM!  14) 

D  WRITE (5. 50) 

D50  FORMAT (  ’  ENTERED  SUBROUTINE’  ) 

BET1 1  .01)  00 
DO  10  1 = I , N 
X( I ) =0 . OD  00 
io  .iriifK  i )  =  i 

0  FIND  THE  ELEMENT  OF  MAXIMUM  ABSOLUTE  VALUE 

DO  I  K  =  I  , N 
C=A!K.K) 

I  1=K 
JJ  =  K 

DO  2  J  =  K ,  N 
DO  2  I  =K,  N 

IF  ( DABS ( C ) -DABS( A( 1 ,J)>  >3,2,2 

3  C  =  A  (  I  .  J  ) 

II  =  1 
JJ  =  J 

2  CONTINUE 

I)ET=DF.T*C 

D  WRITE! 5,55) DET 

055  FORMAT! 5X, *  DET=’D24.17) 

IF! DABS (DET) ) 20, 20, 00 
20  WHITE! 5. 100) 

100  FORMAT!’  MATRIX  A!N,N)IS  SINCULAR* ) 

CALL  EXIT 

30  11!  I  I  >  =  B<  I  I  )/C 

C  DIVIDE  EACH  ELEMENT  OF  THE  IIT1I  ROW  BY  C 

Kf’0  =  K+  I 
DO  4  .1  =K ,  N 

4  A! I  I , J ) = A ( I  I  ,0  )  /C 
A!  1  I  ,  JJ  )  =  1  .01)  00 

IF  (II  . EO.KlCO  TO  60 

0  SWITCH  THE  KT11  now  AND  THE  IITH  ROW 

DO  5  J  =  K  .  N 
C  =  A  ( K  ,  J  ) 

A  ( K ,  .1 )  =  A  ( I  I  ,  J  ) 

5  A  (  I  I  ,  J  )  =  0 
C=R(K) 

»( 10  =  11!  I  l  ) 

B(  I  I  )=C 

60  IF  (JJ  .  KO.KICO  TO  70 

C  STORE  THE  LOCATION  OF  THE  MAX  PIVOT 

I  I  =  Jl’RM!  JJ  ) 

JPRM! .)  J  )  =  Jl’RM!  K) 

Jl’RM  (  K  )  =  1  I 

C  SWITCH  THE  Kill  AND  THE  JJTH  COLUMNS 

DO  6  1  =  1,  N 
C  =  A  <  I  ,  E  ) 

A (  I  .  K  )  =  A ! I  .JJ) 

6  A  (  I  ,  .1.)  )  =  C 

■:  cet  sur  column  of  zeroes  in  the  kth  column 

70  IF<i:.EO.N)CO  10  I 

DO  7  I  =  LI'O  .  N 
DO  0  J  =  F  l’O  ,  N 

O  A!  I  ,  .1  >  =  A  <  1  ,.l  >  -  A  (  I  ,i:>*A(K,J) 

IX  I  )  =  IX  I  )  -  A  «  I  ,  K  U1XK) 

7  A<  I  ,E  ) -»  .  oil  oo 

I  CONTINUE 

C  OBTAIN  IIIE  X’S  IX'  RACK  SUBSTITUTION 

DO  ')  1  =  1  ,  N 
E  N  -  1  +  1 
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24 


Kl’0=K+ 1 

svn=o.o n  oo 

I  F(KPO-N  >24,24,26 
1)0  10  J  =  KI’0,N 
SUN=SUri+X(J>*A<K,J> 

10  CONTINUE 

26  X<K) =B(K) -SUM 

9  CONTINUE 

C  ORDER  THE  X’S 

DO  11  1=1 ,N 

11  BU)  =  X(I) 

DO  12  1=1 ,N 
1 1 = JPRM( 1 ) 

X( 1 1 )  =  B< I > 

12  CONTINUE 

C  WHITE  OUT  THE  VALUES  OF  X(N) 

DO  101  J=1 ,N 

WRITE (5 ,99)J,X(J) 

99  FORMATS, 5X,’  X(  ’  ,  12,  ’  )  =  ’  ,  D24 . 17) 

101  CONTINUE 

RETURN 
ENI) 
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SUBROUTINE  SELECTIN , A, B, VALUE) 

C  SET  UP  CATELOC  ARRAY  ICAT(N) 

DIMENSION  I CAT (14) 

DOUBLE  PRECISION  A (14 . 1 4 )  , B!  14) , VALUE<  14 ) 

C  INITIALIZE  11  ( N )  ,  I  CAT  ( N ) ,  AND  VALUE  !N) 

DO  10  J=  1  ,N 

B( J ) =0 . 0I)+00 
I  CAT ( J  >  =  1 
VALUE (J  >=0. OD+OO 

10  CONTINUE 

C  SELECT  HIE  COEFFI ECl ENTS  TO  BE  INITIALIZED 

WHITE (5, 20) 

20  FORMAT!/, 1  HOW  MANY  COEFFICIENTS  TO  BE  INITIALIZED?’) 

READ! 5 , 30 ) NN 
30  FORMAT! 12) 

DO  60  1=1 , NN 

WRITE! 5, 40) 

40  FORMAT!/ , '  WHICH  COEFFICIENTS  TO  BE  INITIALIZED?’) 

READ! o ,50)11 

50  FORMAT! 12) 

I  CAT! I  I ) =0 
WRITE! 5, 55) 

55  FORMAT!/,*  WHAT  VALUE  DO  YOU  WANT  COEFFICIENTS  TO  BE?’) 

READ! 5 , 5H ) VALUE! II) 

50  FORMAT ( 1)24  .  I?) 

60  CONTINUE 

IN  ORDER  FOR  THE  SET  OF  EQUATIONS  TO  BE  INDEPENDENT  INSERT 
THE  APPROPRIATE  VALUES  IN  BIN)  AND  A(N,N) . 

DO  00  J=l ,N 
DO  70  K= 1 , N 

IF! ICAT(K) .NE.0)GO  TO  70 
A(K, J)=0. OD+OO 
A!K.K)=1 . 0D+00 
BIK) =VALUE!K) 

70  CONTINUE 

00  CONTINUE 

DO  90  J=1  ,N 
WRITE! 5 , 1 040 ) J , B! J) 

1040  FORMAT!/, ’  DC ,12,’ ) = ’ , D24 .17) 

90  CONTINUE 

RETURN 
END 
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SUBROUTINE  SUI'RES(N,A,B,X> 

C  SET  UF  CATELOC  ARRAY  JCAT(N) 

DIMENSION  JCAT (14) 

DOUBLE  I'UECISION  A(  14 . 14 )  ,B<  14)  ,X(  14) 

C  INITIALIZE  JCAT(N)  TO  ONE 

DO  20  J  =  1  ,  N 

J  CAT ( J )  =  I 

20  CONTINUE 

WRITE!  5,40) 

40  FORMAT!/, ’  HOW  MANY  VALUES  OF  JCAT=0’ ) 

HEAD ( 5 , 00 )KK 
60  FORMAT (12) 

DO  00  L= l ,KK 

WRITE (5, 90) 

90  FORMAT!/, 1  WHICH  VALUES  OF  JCAT=0* ) 

READ (5, 95) I 

95  FORMAT  (12) 

JCAT< I ) =0 

00  CONTINUE 

C  SET  TI1E  ROW  IN  MATRIX  A ( N , N )  =  0 

DO  300  J=1 , N 
DO  200  K=  1  ,  N 

IF(JCAT(K) .NE.0JGO  TO  200 
150  A(K,J)=O.0D+0O 

A(K,K)=1 .0D+00 
BOO  =0 .0D+00 

D  WRITE(  5 , 1  00  )K ,  .1 ,  A(K,  J  ) 

I)  1 00  FORMAT!/’  AO, 12,’  ,,,I2,’  )=’,D24.17) 

200  CONTINUE 

:>()0  CONTINUE 

RETURN 
END 
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D 

1)1095 

C 


1) 

1)1  ICO 
20 

21 

24 


31 

35 

C 


CO 

I) 

1)1  125 
C 


0 


c 


r.o 

1030 


SUBROUTINE  YCALCU !N,X, YCALC, XI NC) 

DOUBLE  PRECISION  X(14) 

DIMENSION  XI  NCI  500) ,YCALC(500> .XSING! 14) ,X( 1 ) ,Y(1 ) 

DATA  DENI , BEN2 , DEN3 , DEN4 , DENS ,DEN6 ,XNUMl ,XNUM2  ,XNUM3  ,XNUM4 
1 XNUM5 , XNUM6/ 1 2*0 . 0/ 

BYTE  FI I. EX (32) 

VRITEC5, 1095) 

FORMAT ( 1  ENTERED  SUBROUTINE  YCAL’ ) 

PUT  X(N)  INTO  SINCLE  PRECISION 
DO  20  L=  1  ,  N 

XS I NO! L) =SNCL! X( L) ) 

WRITE ( 5 , 1 IOO)XSING(L) 

FORMAT!/’  XS  I  NC= '  ,  El  5 . 5 ) 

CONTINUE 
WRITE! 5 , 2 1 ) 

FORMAT!/ , ’  MAXIMUM  X  YOU  WISH  TO  CO’ ) 

READ (5,24) IX 
FORMAT ( 14) 

J= IX+l 
WRITE! 5 , 3 i ) 

FORMAT ( / , '  X  INCREMENT  VALUE?’ ) 

READ ( 5 , 35 )XI WCRE 
FORMAT ( F 1 0 . 0 ) 

ALLOCATE  STORAGE 
DO  30  I  =  1  .  J 

XINC( I )-0.0 
YCALC( I )=0.0 

CONTINUE 
WRITE! 5, 1 125) 

FORMAT ( ’  HAVE  SET  UP  STORACE’ ) 

SET  UP  LOOP  TO  CALCULATE  YCALC(K) 

DXI NC=0 . 0 
DO  50  K=  1  ,  J 

X1NC(K)=DXINC 
DX I NC2  =  DX I NC*DX I NC 
I)X  1  NC3  =  DXI  NC2*DXI  NC 
DX I NC4=  DX I NC2*DXINC2 
DXI  NC5  =  I)XI  NC2*I)XINC3 
I)X  I NC6  =  DX  I  NC3*DX  I NC3 
DEN  I =  XS I NG ( 2 ) *DX I NC 
D£N2=XSINC(3)*DXINC2 
DEN3=XS 1  NO! 4 ) *DXI NC3 
DEN4  =  XS I NGC  5 )*DXI NC4 
DEN5=XSING(6)*DXINC5 
DEN!»  =  XS1  NC(7  )*DXINC6 
XNUM1 =XS I NC< 9  >*DXI  NC 
XNDH2=XSINCC I0)*DXINC2 
XNUM3=XS I NC( t 1 )*DXINC3 
XHUM4=XSINC( I2)*DXINC4 
XriUM5=XSIHC(  1 3)*I)XINC5 
XN1JM(»=XSING(  1 4  )*DXI  NC!> 

FORM  NUMERATOR  AND  DENOMINATOR  OF  YCALC 

XN1JM  =  XS  I  NC((J)  +XNUMI  +XNUH2+XNUM3+XNUM4+XNUM5+XNUM6 
XDEN  =  XS I NC( I )+I)EN1 +DEN2+DEN3+DEN4+DEN5+DEN6 
CALCULATE  YCALC 

YGAI.C(K)  =XNUM/XI)EN 
INCREMENT  DXINC 

DXI NC=DXI NC+X1NCRE 

CONTINUE 
WRITE! 5. 1050) 

FORMAT!/’  *F1 RENAME  FOR  X.Y  DATA:  ’) 

REA D ( 5 , 1 060 )  I  .EN , F I  LEX 
FORMAT (0.32A I ) 

F I  LEX ! LEN  +  1)*0 


1 050 


0PEN(UN1T=!  ,NANE=  FI  LEX,  FORM=*  UNFORMATTED’  ,TYPE=’1TEW’  ) 
DO  220  1=1 ,J 
READ< 1)  X( l ) ,Y( 1 ) 

220  CONTINUE 

CI.OSE(UNIT=  1  ) 

CONTINUE 

RETURN 

END 
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APPENDIX  E 


DETERMINATION  OF  THE  VOLUME  OF  A 
SAMPLE  BY  ARCHIMEDES'  PRINCIPLE 


An  accurate  method  to  determine  the  volume  of  a  sample  is  by  the 
use  of  Archimedes'  principle.  The  volume  of  a  sample  at  ambient 
pressure  P^  and  room  temperature  Tq  (%25°C)  is  given  by 


V(P  T  )  = 
o  o 


*)(*  -  &) 


H20  -  Pa 


To  calculate  the  volume  of  a  sample  using  Eq.  (El),  one  must  determine 

the  mass  of  the  sample  mg  in  air,  determine  the  mass  of  the  sample  In 

water  m^  Q,  and  use  a  known  value  of  the  density  ph  Q  of  the  water, 

the  density  of  the  weights  used  on  the  balance,  and  the  density 

of  the  air  in  the  room  at  the  temperature  T  in  the  room.  The  masses 

o 

in  Eq.  (El)  are  apparent  masses. 


In  this  investigation,  the  masses  m^  and  Q  were  determined  by 

using  a  Mettler  Model  H-311  Triple-Beam  Balance,  which  has  a  maximum 

capacity  of  240  g  and  is  accurate  to  ±0.0003  g.  Distilled  water  was 

used  in  this  volume  determination.  It  was  deaerated,  covered,  and 

then  cooled  to  25°C.  The  density  P„  _  of  the  distilled  water  [17] 

n^U 

z  3 

and  the  density  p  of  the  air  were,  respectively,  0.99705  g/cm  and 

3 

0.001185  g/cm  .  An  average  from  four  volume  determinations  was  taken 
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for  both  the  butyl-252  sample  and  the  type-W  neoprene  sample.  Their 

3  3 

volumes  were,  respectively,  8.8433  cm  and  8.7774  cm  . 
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APPENDIX  F 

A  FORTRAN  PROGRAM  (PYCNOM.FTN)  TO  CALCULATE  THE 
SAMPLE  VOLUME  AS  A  FUNCTION  OF 
TEMPERATURE  AND  PRESSURE 


The  FORTRAN  program  PYCNOM.FTN  calculates  the  volume  of  the  sample 
used  in  the  acoustic  pycnometer  according  to  Eq.  (14b).  The  listing 
of  the  program  is  given  as  follows: 
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rtnftnno 


1100 
1 1 10 

195 

190 

194 

198 

199 
190 

105 

1185 

1200 

197 

1199 

1201 

070 

C 


10 

G 


20 

G 

25 
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THIS  PROGRAM  CALCULATES  THE  VOLUME  AND  THE  ISOTHERMAL 
BULK  MODULUS  AS  A  FUNCTION  OF  PRESSURE. 

THIS  PROGRAM  STORES  THE  PRESS.  AND  VOLUME  IN  FILE 
CALLED  * FOR0O2 . DAT’ 

BYTE  ANS14) 

DIMENSION  A(14).AC!14),C(6) .PRESS! 100) .VSAMPI 100) 

1 , DERI  V ( 100 ) .BULK! 100) 

WHITE! 5 , 1 100) 

FORMAT!/.’  SUSE  OLD  DATA?’) 

READ! 5 , 1 1 10.END=8Q0)ANS 
FORMAT (4A1 ) 

IF  (ANSI  1 ) .EQ>. ’ N’ )CO  TO  870 

OPEN(UNlT= 1 , NAME= ’ PYCNOM.DAT’ ,TYPE=’OLD’ ) 

READ(1 , 190) VO 
WHITE! 5 , 1 95 ) VO 
FORMAT!/,’  V0= ' ,F!2. 5) 

READ!  1  , 190)TEM1’ 

WRITE! 5 , 1 96 ) TEMP 
FORMAT!/.’  TEMP  =’,F12.5) 

READ!  1  ,  1 90  )IU!0 
WRITE!  5. 194)11110 
FORMAT!/,  ’  R1I0=  ’  , F12 . 5) 

READ! 1 . 190)WTSR 
WRITE! 5. 198)WTSR 

FORMAT!/.’  WEIGHT  OF  I1C  FOR  SAMPLE  RUN  =’ .F12.5) 

READ! I , I90)WTCR 
WRITE! 5, 199)WTGR 

FORMAT!/,  ’  WEIGHT  DURING  CALI B.  =’  .F12.5) 

FORMAT! El 2 . 6 ) 

READ! 1 ,  1G5)NUMCR 
FORMAT! 12) 

DO  1200  K= 1 , NUMCR 
READ!! , 1 185 ) AC(K) 

FORMAT! El 2 . 5 ) 

CONTINUE 

READ! 1 , 1 97 )NUMSR 

FORMAT! 12) 

DO  1201  J  =  l  , NUNSR 
READ! 1 , l 199)A(J) 

FORMAT! El 2. 5) 

CONTINUE 
GO  TO  1210 

OPF.N!UNIT=  1  ,NAME=’ PYCNOM.DAT’  ,TYPE=’NEW’  ) 

ALLOCATE  STORAGE 
DO  10  J=t , 100 

PRESS! J) =0.0 
VSAMP ! J ) =0 .0 
DERIV! J >=0.0 
BULK!.) )  =0 .0 

CONTINUE 

DO  20  1=1,14 

A! ’  )=0.0 
AC! I ) =0 .0 

CONTINUE 

DO  25  K=l  ,6 

C !  K )  =  0 . 0 

CONTINUE 
WRITE (5. 30) 

FORMAT!/.  ’  CENTER  THE  I NITI AL  SAMPLE  VOLUME’ ) 
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READ(5 , 00  >  VO 
WRITE! 1 , 290 ) VO 

290  F0RMAT(F12 . 5 ) 

VUITE( 5  35) 

35  FORMAT!/,  ’  SENTER  TEMP.  OF  SAMPLE  RUN’ ) 

HEAD ( 5 , 90 )TEMP 
WniTE! 1 ,291 )TENP 

291  FORMAT! F 1 2.5) 

WRITE! 5 , 36) 

38  FORMAT!/,’  SENTER  DENSITY  OF  HG’ ) 

READ!  5 , 90  )R1I0 
WRITE!  I  ,  292)11110 

292  FORMAT! F 12. 5 > 

WRITE! 5 , 40 ) 

40  FORMAT!/, ’  SENTER  WT.  OF  HG  DURING  SAMPLE  RUN’) 

READ! 5 , 90 )WTSR 
WRITE! 1 , 293 )WTSR 

293  FORMAT! FI  2. 5) 

WRITE  (5,50) 

30  FORMAT (/. ’  SENTER  WT.  OF  HG  DURING  CALIBRATION’) 

READ! 5 , 90 )WTCR 
WRITE! I , 294)WTCR 

294  FORMAT ( F 1 2 . 5 ) 

90  FORMAT! El 2. 6) 

C  GET  READY  TO  READ  THE  COEFFICIENTS 

WRITE! 5,60) 

60  FORMAT!/. ’  SENTER  THE  NUMBER  OF  COEFFICIENTS  IN  CALIB.  FORMULA’) 
REA  D  ( 5 . 6  5 )  NUf ICR 
WRITE! 1 ,65)NUMCR 

65  FORMAT! 12) 

DO  95  J  = 1 , NUMCR 

WRITE! 5 ,80)J 

80  FORMAT!/,’  SENTER  AC! ’, 12, ’)’ ) 

READ (5,85 ) AC ( J ) 

WRITE! I , 85 ) AC! J ) 

05  FORMAT! El  2. 5) 

95  CONTINUE 
WRITE! 5 ,96 ) 

96  FORMAT!/, ’  SENTEIl  THE  NUMBER  OF  COEFFICIENTS  IN  THE  SAMPLE  RUN*) 
READ! 5,97) HUMSR 

WRITE! 1 , 97 )NUMSR 

97  FORMAT! 12) 

DO  90  K=  1  , NUMSR 

WRITE! 5, 99 )K 

99  FORMAT!/,’  SENTER  A! ’ , 12, ’ ) ’ ) 

READ! 5, 101 )A(K) 

WRITE! 1 ,101 ) A(K) 

101  FORMAT ( El  2 . 5 ) 

98  CONTINUE 

1210  CLOSEOIN IT=  I  ) 

0  INITIALIZE  VARIOUS  VALUES  TO  ZERO 

VSNUM=0 ,0 
DVSNI)M=0 . 0 
VS DEN =0 . 0 
DVSREN-O.O 
VltSIUO.O 
TERM1 =0.0 
TEIIM2  =  0 . 0 
DI FF=0 .0 
CDEN=0 .0 
DERI  V 1=0.0 
DERI V2= 0.0 
VCNUN=0.0 
DVCN1IM=0.0 


VC1)KN=0 . 0 
DVC1)EN=0.0 
TERMA=0 .0 
TKRMU=0 .0 
FACTOR® 0 . 0 
FACT 1=0.0 
1>EUFAC=0.0 
Vl*T=0 . 0 
DERIV3=0.0 

vucn=o.o 

P  =  0.0 
P2=0.0 

pa=o.o 

P4=0 . 0 
P5=0.0 
P6=0 .0 
P 1 HC=0 .0 

ENTER  THE  VALUES  OF  C<6> 

C(I )=0. 1000E+01 
C!2)=0. IG21E-03 
C!3)  =  -0.3753E-!)4 
C  <  4 )  =  -0 . 5 1 9E-07 
C!5)  =  0 .  454E-08 
C!6)  =  0. 1 1 9E- 1 0 

START  PERFORMING  THE  CALCULATIONS 
TO=25 . 0 

CDEN=C( 1 )+C(2)*T0 
DIFF=!WTCR-WTSR)/RHO  -  VO 
FACTOR® D I FF/CDEN 
WRITE! 5, 120) 

FORMAT ( / , ’  3 ENTER  MAXIMUM  PRESSURE’) 

REAIM  5, <U0)KP 
FORMAT! 14) 

WRITE! 5, 140) 

FORMAT!/,’  $ ENTER  PRESSURE  INCREMENT’) 

READ! 5 , I  30 ) PI NC 
FORMAT! FI 0.0) 

INITIALIZE  THE  PRESSURE  TO  ZERO 

P=0.0 

JP=KP+ 1 

ALLOCATE  STORACE 
DO  1500  J=I  ,  JP 

PRESS! J) =0.0 
VSAMP! J ) =0 . 0 
DERIV! J)=0.0 
BULK! J )  =0 . 0 

CONTINUE 
DO  100  l=l,JP 

PRESS! I )=P 

CALCULATE  P2 , P3 , P4 , P5 , P6 

P2=P*P 

P3=P2*P 

P4=P2*P2 

P5=P3*I*2 

P6  =  I’3*P3 

CALCULATE  TIIE  VARIOUS  VOLUMES 

CALCULATE  VOLUME  OF  BORE  FOR  THE  SAMPLE  RUN 

VSNUH=A!a)+A!9)*P+A< 10)*P2+A! 1 1 >*P3+A< I2)*P4+ 

I  A!  13)*PS4A(  I4)*P6 

CALCULATE  THE  DERIVATIVE  OF  V-'NUH 

DVSHUM® A <  9 ) +2 . 0*A 1 1 0 ) *P+3 . 0*A< 1 1 ) *P2+4 . 0*A ( 1 2 )*P3 
+5.0*A! I3)*P4+6.0*A! I4)*P5 
VSDEN  =  A!  1  )+A<2  )*l’+A!  3  )^P2+A!4  )*P3+A(  5)*P4-*-A!6)*P5 
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1  +A(7)*P6 

CALCULATE  THE  DERIVATIVE  OF  VSDEN 

DVSDEH  =  A ( 2 ) +2 . 0*A ( 3 ) *P+3 . 0*A ( 4 ) *P2+4 . 0*A ( 5 ) *P3+ 

I  5 . 0*A<  6 )*P4+6 . 0*A(7 ) *P5 

CALCULATE  THE  VOLUME  OF  THE  BORE  DURING  TUE  SAMPLE  RUN 
VBSR=VSNUM/VSDEN 

TERM1=VBSR/VSNUM 
TEllf  12  =  V  I!SR/  VSDEN 

CALCULATE  THE  DERIVATIVE  OF  VBSR 

DEll  I  V  1  =TEHM!  *DVSNUM-TERM2*DVSDEN 
START  CALCULATIONS  FOR  CALIBRATION  RUN 

VCNUM= AC ( B ) +AC ( 9 ) *P+AC (10) *P2+AC (11) *P3+AC ( 1 2 )*P4 
1  +AC( I3)*P5+AC( I4)*P6 

DVCNUM=AC(9 ) +2 . 0*AC( I  0 ) *P+3 . 0*AC' 1 1 )*P2+4 . 0*AC< 12)*P3+ 

1  5 . 0*AC( 1 3 )*P4+6 . 0*AC( 1 4 )*P5 

VC1)EN  =  AC  ( I  )  +  AC  ( 2 )  *P+ AC  ( 3 )  *P2+AC  ( 4 )  *P3+AC  ( 5 )  *P4+AC  ( 6  )*P5 
1  +AC(7)*P6 

DVCDEN  =  AC ( 2 ) +2 . 0*AC ( 3 ) *P+3 . 0*AC ( 4 ) *P2+4 . 0*AC ( 5 ) *P3+ 

1  5 . 0*AC ( 6 ) *P4+6 . 0*AC(7 )#P5 

CALCULATE  THE  VOLUME  IN  THE  BORE  FOR  THE  CALIBRATION 
VBCR= VCNUN/VCDEN 

TERMA=VBCR/VCNUM 
TERM U=V BCR/ VCDEN 

CALCULATE  THE  CORRESPONDING  DERIVATIVE 

DERI  V2=TERHA*DVCNUM-T) RMB*DVCDEN 
CALCULATE  THE  VOLUME  DUE  TO  A  DIFFERENT  FILL  HEIGHT 

FACT1=C( 1 ) +C(2 )*TEMP+C(3 )*P+C(4 )*TEMP*P+C( 5 )*P2+ 

1  C  ( 6 )*TEHP*P2 

V  PT  =  F  ACT 1  * F ACTOR 
CALCULATE  THE  DERIVATIVE  OF  VPT 

I)ERFAC=C( 3 ) +G( 4 ) *TEMP+2 . 0*C ( 3 ) *P+2 . 0*C ( 6 )*TEMP*P 
DERI V3=DERFAC* FACTOR 

CALCULATE  THE  SAMPLE  VOLUME  AND  IT’S  DERIVATIVE 
VSAMP ( I )=VO+VBSR-VBCR+VPT 
DER I V ( I ) =  DERI V I -DERI V2+DERI V3 
CALCULATE  THE  COMPRESSIBILITY 

CONI’*  -  ( I  .  0/VSANP(  I  ) )  *  DER  I V  ( I  ) 

CALCULATE  THE  ISOTHERMAL  STATIC  BULK  MODULUS 
BULK ( I )= I .0/COMP 
INCREMENT  THE  PRESSURE 
P=P+PINC 

CONTINUE 

WRITE  OUT  THE  VOLUME  AS  A  FUNCTION  OF  PRESSURE 
DO  200  .1  =  1  ,JP 

WRITE ( 5 , 1 60) PRESS (J ) , VSAMP (J ) 

WRITE (2 , 9999 ) PRESS ( J ) , VSAMP (J ) 

F0RMATI2E12.5) 

FORMAT</.'  PRESSURE  (MPA)  =  ’ . El 2 . 5 . 3X. ’ VSAMP  (ML)*’ ,E12.5> 

WRITE  OUT  THE  BULK  MODULUS 

WRITE! 5, 170 ) PRESS (J ) ,BULK(J) 

FORMAT!/ , ’  PRESSURE  ( MPA )  =  ’  , E 1 2 . 5 , 3X , ’ BULK  MOD .  (MPA)  =  •  ,EI2. 5 ) 

CONTINUE 

CI.0SE(UHIT=2) 

WRITE! 5. 1000) 

FORMAT ( / , *  PROGRAM  EXITS!’) 

CALL  EXIT 
END 
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APPENDIX  G 


PICTURE  OF  ACOUSTIC  PYCNOMETER  AND  ENGINEERING  DRAWINGS 


Fig.  G1  —  Assembled  acoustic  pycnometer 
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#5  (.205)  DRILL  THRU 
C'BORE  5/16  DIA.  x 
1/4  DEEP 


Fig.  G2  —  Thermistor  plug 


APPENDIX  H 


STATIC  BULK  MODULUS  CALCULATION 

The  static  bulk  modulus  of  the  elastomer  samples  used  in  the- 
acoustic  pycnometer  measurement  can  be  calculated  from  the  least- 
squares  coefficients  obtained  by  fitting  the  experimental  data  to  a 
sixth-degree  polynomial.  The  sixth-degree  polynomial  expression  for 
the  volume  of  the  elastomer  sample  as  a  function  of  pressure  at  constant 
temperature  is  given  by 


V SAMPLE  '  A!  +  *2P  +  V2  +  V’  +  +  V’  +  *7p6  '»*  > 

where 

Aj  through  A^  -  the  least-squares  coefficients  of  the  sixth-degree 
polynomial 

The  Isothermal  compressibility  k  Is  obtained  by  taking  the  derivative 
of  Eq.  (HI)  with  respect  to  pressure  and  dividing  this  by  the  volume  of 
the  sample  at  the  desired  pressure.  The  Isothermal  compressibility  is 


(H2) 


The  static  bulk  modulus  B  is  the  reciprocal  of  the  isothermal 
compressibility  *.  The  expression  used  to  calculate  the  static  bulk 
modulus  B  at  constant  temperature  is 
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i  pi“ 

B(P)  -  -  — ~  -  .  (H3) 

y^a  -  oa^1-2 

i  =  2 

The  static  bulk  modulus  of  the  butyl-252  sample  at  10°C  and  25°C 
is  shown  in  Fig.  HI.  (The  behavior  of  the  curves  above  65  MPa  is 
believed  to  have  no  physical  significance.  This  behavior  seems  to 
be  an  artifact  of  the  least-square  analysis) . 


98 


REFERENCES 

1.  R.J.  Roark  and  W.C.  Young,  Formulas  for  Stress  and  Strain, 

New  York:  McGraw-Hill  Book  Co.,  5th  Edition,  pp.  4,  1975. 

2.  R.J.  Bobber,  Underwater  Electroacoustic  Measurements,  Washington, 
D.C.:  U.S.  Government  Printing  Office,  pp.  312,  1970. 

3.  A.J.  Rudgers,  "Theory  of  a  Method  for  Measuring  the  Complex  Dynamic 
Bulk  Modulus  of  Elastomers  by  use  of  Rayleigh-Scattering,"  Draft 

of  NRL  Formal  Report  to  be  published. 

4.  J.W.S.  Rayleigh,  The  Theory  of  Sound,  New  York:  Dover  Publications 
Vol.  II,  pp.  283,  1945. 

5.  R.D.  Corsaro,  J.  Jarzynski  and  C.M.  Davis,  Jr.,  "Acoustic 
Densitometer  with  Results  for  Polyethylene  Oxide,"  Journal  of 
Applied  Physics,  Vol.  45,  No.  1,  pp.  1-8,  1974. 

6.  Ibid 

7.  Ibid 

8.  "CRC  Handbook  of  Chemistry  and  Physics,"  R.C.  Weast,  ed..  West 
Palm  Beach,  FL:  CRC  Press,  Inc.,  59th  Edition,  pp.  F-ll,  1978-1979 

9.  E.G.  Kimmich,  "General  Engineering  Properties  of  Rubber,"  Chap.  4, 
Engineering  Uses  of  Rubber,  A.T.  McPherson  and  A.  Klemm,  eds.. 

New  York:  Rheinhold  Publishing  Corporation,  pp.  72-100,  1956. 

10.  Ibid 


II.  Matec,  "Special  Instruction  Manual  for  Using  Model  122-B  in 
Multiple  Instrument  System,"  Matec,  Inc..  Warwick,  R1 
i!.  ibid 


1U1 


«.v„ 


REFERENCES  NOT  CITED 

V.C.  Anderson,  "Sound  Scattering  from  a  Fluid  Sphere,"  The  Journal 
of  the  Acoustical  Society  of  America,  Vol.  22,  #4,  pp.  426,  July  1950. 
J.J.  Faran,  "Sound  Scattering  by  Solid  Cylinders  and  Spheres," 

The  Journal  of  the  Acoustical  Society  of  America,  Vol.  23,  #4,  pp.  405, 
July  1951. 

P.M.  Morse,  Vibration  and  Sound,  New  York:  McGraw-Hill  Book  Company, 
Inc.,  1936. 

S.N.  Rschevkin,  The  Theory  of  Sound,  New  York:  The  Macmillan  Co., 
Pergamon  Press,  Chapter  9,  1963. 

R.W.  Timme,  "Sound  Speed  in  Castor  Oil,"  The  Journal  of  the  Acoustical 
Society  of  America,  Vol.  52,  #3,  pp.  989,  1972. 


103 


